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3.0 ANALYSIS OF INSTRUMENTAL POLARIZATION OF IMAGING OPTICS 
3.1 Introduction 

The objective of the SAMEX magnetograph's optical 
system is to accurately measure the polarization state of 
sunlight in a narrow spectral bandwidth over the field of view of 
an active region to make an accurate determination of the 
magnetic fields in that region. Our design goal is to measure 
magnetic fields to an accuracy of one part in 10 4 . To achieve 
this accuracy requires a polarimetric accuracy of 10“ 4 in 
determining the polarization components of the light as a 
fraction of the total intensity. This requirement means that the 
instrumental polarization of the optics must be reduced to levels 
below 10~ 5 . 

All optical elements introduce some polarization 
change/ especially when used off axis. Combinations of mirrors 
and antireflection-coated lenses can display a full range of 
polarization behavior: linear and circular polarization and 

linear and circular retardancs. Qf particular concern in the 
development of highly accurate polarimeters is polarization 
rotation which causes linearly polarized light to leak through 
subsequent of crossed polarizers. 

Because of this instrumental polarization, the SAMEX 
foreoptics (the optical elements in front of the polarizer - the 
Cassegrain telescope and the relay lenses) must be considered as 
a weak polarizer in front of the polarimeter section. This 
instrumental polarization in the foreoptics changes the 
polarization state of the sunlight incident on the polarimeter 
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and thus introduces errors in the measurement of the. solar 
maanetic field. 

If possible, it would be preferable to place the 
polarimeter before the imaqina optics, to locate polarizers and 
retarders in front of the Casseqrain telescope. Then the 
determination of the polarization state of liqht would be 
unaffected by the instrumental polarization of the foreootics. 

t 

This desiqn would have the additional advantaqe of allowinq the 
liqht to pass throuqh the polarimeter over a smaller ranqe of 
anqles of incidence. However, such a desiqn is impractical for a 
system with a 30 cm aperture - hiqh quality polarizers and 
retarders have much smaller apertures, on the order of a few 
centimeters. Thus, it is necessary to use collect inq optics to 
collect the 30 cm aperture of liqht and focus it through small 
polarizing elements in the polarimeter section. As the collected 
liqht is passed throuqh smaller apertures, the anqular spectrum 
of the liqht increases in a relationship governed by the Lagrange 
invariant. By reducing the beam from 30 cm to 2 cm throuqh the 
tunable filter and polarimeter, the range of anqles of incidence 
is increased by a factor of 15, from 5 min of arc to 75 min of 
arc, thereby increasing the anqle of incidence effects in the 
lenses just in front of the polarimeter and in the polarimeter 
itself. Consequently we anticipate there will be some 
instrumental polarization due to the foreoptics in the liqht 
entering the polarimeter. 

The goal of the analyses outlined in this section is to 
precisely characterize the extent of this instrumental 
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polarization, and to design the optics and coatings in unison to 
minimize this spurious polarization introduced by the 
foreoptics. We will show calculations of the instrumental 
polarization of ordinary foreoptics compared with the 
polarization performance of special ultra-low polarization 
optical coatinas designed for this application. 

The instrumental polarization analysis uses a program 

* 

tha incorporates the theory of polarization into the standard 
geometrical optics and lens design codes (Chipman, 1987). By 
including the polarization of the optical elements in the first- 
and third-order design process, the effects of coatings on curved 
substrates can be treated. For each ray, a polarization matrix 
(in the Jones matrix formulation) is calculated for the ray at 
each optical interface. These matrices are multiplied together 
to calculate the polarization matrix for that ray from object 
space to image space. An analytic function for the rays provides 
the polarization behavior as a function of the exit pupil, object 
height, and wavelength. This technigue represents a quantum jump 
improvement in the practical desiqn of foreoptics in front of a 
polarimeter and will enable the degree of residual instrumental 
polarization to be reduced to < 10“ 5 in the SAMEX magnetograph 
system. 


3.2 Instrumental Polarization of Standard Cassegrain Telescopes 
To establish the need for a detailed polarization 
analysis, we estimate the polarization effects associated with a 
Cassegrain telescope with aluminum thin film coatings. The 
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electric vector of light incident on a surface of the Cassegrain 
can be decomposed into two parts, the components of the vector 
vibrating parallel (p) and perpendicular (s) to the plane of 
incidence. These components have different reflections as a 
function of the angles of incidence. Figure 23. a shows the 
reflectance for s and p polarized light from an aluminum thin 
film coating with complex index of refraction n = 0.7 - 7.0 T. 
Figure 23. b shows the phase change on reflection for the s and p 
components. The difference in s and p reflectance causes a weak 
linear polarization aligned perpendicular to the plane of 
incidence to be associated with reflection from all metallic 
interfaces. Moreover, the differences in phase change cause a 
weak linear retardance to be associated with reflection from 
mirrors . 

The percent reflections from the aluminum surface for 
the surface parallel (s) and surface perpendicular (p) rays are 
approximately given by 

R s = 0.94 + 0 . 02i 2 , 


and 


R p = 0.94 - 0 . 0 2 i 2 , 


where the reflectance rate of change per 
squared (i 2 ) is approximated from Figure 
induced linear polarization is estimated 


unit angle of incidence 
23. a. Therefore, the 
to be 
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Figure 23. Linear polarization and retardance effects for a 
Cassegrain telescope. The percent reflectance (a) and phase 
change (b) on reflection from an aluminum thin-film coating is 
shown for both the s and p components of the incident light wave. 
The wavelength of the light is 5250 A and the complex refractive 
index of the aluminum coating is n * 0.7 - 7.0 T. The 
differences in the s and p reflectance cause linear polarization 
aligned with the incident plane. The differences in the s and p 
phase cause a linear retardance. The differences are small for 
.small angles of incidence, but they are. not negligible for the 
^.'SAMEX magnetograph design. 




0 . 02i 2 . 



At the edge of the primary mirror of aperture D and focal length 
£ 1 , a ray is reflected through a total angle of tan(D/fl) ~ D/fl 
- 1/f, where f is the f-ratio. Then for a Cassegrain telescope 
the maximum angle of incidence i m is given approximately by one 
half the reciprocal of the f-ratio. For an f-ratio of 5, the 
angle of incidence is i m *l/10. Hence the induced linear 
polarization is on the order of 3 2 x 1Q~ 4 . For a linearly 

polarized ray this represents a rotation of its plane of 
polarization over a segment of the mirror, and this rotation 
introduces errors in the deduced polarization state. In Figure 
24 the angle of incidence versus the pupil coordinate for an 
'illustrative' Cassegrain telescope is shown for an on-axis and 
off-axis ray. The point to note is that, for the off-axis ray, 
the average angle of incidence is not zero, and hence there is a 
net polarization effect associated with the off-axis rays. 
Furthermore, even on-axis, we average over the square of the 
incident angle which gives a net polarization contribution, as 
shall be discussed below. 

For an on-axis beam incident on either of the aluminum 
coated Cassegrain telescope mirrors, the magnitude and 
orientation of the linear polarization and linear retardance 
associated with the mirrors have the forms shown in Figure 25. a 
(polarization) and Figure 25. b (retardance) as a function of the 
pupil coordinates. In Figure 25. a, the linear polarization is 
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Figure 24. Illustrative example of angles of incidence at 
Cassegrain mirrors. These curves show the variations of the angle 
of incidence along the normalized pupil coordinate for on-axis 
and off-axis rays incident on both the primary and secondary 
mirrors of a Cassegrain system. The off-axis rays have a non- 
*zero average angle of incidence at both the primary and secondary 
mirrors. This implies a net linear polarization associated with 
the telescope for off-axis rays. The example shown here is for a 
Cassegrain with a larger field of view than the one chosen in the 
SAMEX design. It was chosen to provide a vivid example of the 
off-axis problem. 






/ 

\ 

/ 


1 • • 

• 1 

N 

✓ 

\ 

/ 


B. 




% 



_• • . • -.V. - . • c 

TIgure”25T '.Patterns of linear polarization and retardance for -a 
Cassegrain telescope. The linear polarization pattern (a) and' 
linear retardance pattern (b) associated with an aluminum-coated 
Cassegrain primary mirror are shown for an on-axis beam. In (a) 
the magnitude and orientation of the linear polarization are 
given as linear segments. In (b) the linear segments represent 
the magnitude and fast-axis orientation of the retardance. Both 
patterns are quadratic functions of the pupil coordinate and are 
termed "polarization defocus." 





zero in the center of the mirror where the beam is at normal 
incidence. The magnitude of the linear polarization 
(polarizance ) increases quadratically with distance from the 
center of the mirror. The linear polarization is oriented 
tangentially. This polarization pattern (polarization 
aberration) is called linear polarization defocus (or quadpol). 
Figures 26. a and 26. b show the effect of large amounts of linear 
polarization defocus (about 40 times more than the conventional 
Cassegrain) on beams of uniform lef t-circularly polarized light 
(a) and uniform, vertical-linearly polarized light (b). The 
linear polarization associated with the telescope mirrors changes 
the polarization state of the light causing spatial variations of 
intensity and polarization across the beam. If the polarization 
state of the light is now measured with a polarimeter, a 
polarization state different from that incident on the mirrors is 
obtained. Despite the symmetry associated with the resulting 
transmitted polarization patterns (such as in Figure 26), the 
polarization variations do not cancel (due to averaging over the 
incident angle squared). This is best understood by considering 
the transmitted light as being a superposition of two 
polarization states: the incident state yields the correct 
polarization measurement; the light in the orthogonal state 
constitutes the error signal introduced by the instrumental 
polarization . 

.For example, in the case of Figure 26. b, the light in 
the orthogonal state (horizontal linear polarization) has the 
form across the pupil shown in Figure 27. Although the phases of 
this horizontal component are 180° out of phase in the four 
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Figure 26. Linear polarization defocus effects for a Cassegrain 
telescope. This figure illustrates linear polarization defocus 
effects for a Cassegrain primary mirror for a beam of (a) uniform 
left circular polarization and of (b) uniform vertical linear 
polarization. The ellipses and arrows at the top, bottom, left, 
and right of the diagram represent the same relative position of 
the polarization ellipse of the reflected beam of a Cassegrain 
mirror with an exceptionally large field of view (40X the normal 
Cassegrain field). The location of the arrowheads represents the 
1 phase of the light where one cycle is a full wave. The defocus 
, effect introduces the orthogonal state of polarization which 
'represents an error signal introduced by the instrumental 
polarization in the measurements made with the magnetograph. 





f! 





Figure 27. Polarization errors associated with a Cassegrain 
telescope. This figure illustrates the orthogonal state of 
polarization introduced by the instrumental polarization of a 
Cassegrain mirror for the incident linear (vertical) polarization 
of Figure 26 b. These horizontal components will pass through 
the linear polarizer of a magnetograph's polarimeter when it is 
in the horizontal position and thereby contribute to the error 
/signal of the magnetograph. 


quadrants (resembling astigmatism), all this light will pass 
through a horizontal linear polarizer and contribute to the 
polarization error signal. The phase differences do not cause 
cancellation of the polarization aberration; instead they affect 
the polarization accuracy and change the structure of the 
diffraction pattern (Kubota and Inoue, 1959). 

For a metallic mirror, the effects of the linear 
retardance are orders of magnitude larger than the effects of the 
weak linear polarization. Figure 25. b shows the form of the 
linear retardance associated with an on-axis beam incident at a 
Cassegrain telescope mirror. The orientations of the lines 
represent the orientation of the fast axis of the retardance. 

The lengths of the lines signify the magnitude of linear 
retardance, which increases quadrat ically with pupil 
coordinate. This polarization aberration is called linear 
retardance defocus or quadtard. Figures 28. a and 28. b show the 
effect of quadtard on a uniform lef t-circularly polarized beam 
(a) and on a uniform, vert ical— linearly polarized beam (b). The 
dominant effect of the retardance is a coupling of linear into 
circular polarized light and vice versa. These figures are for 
retardances about 40 times greater than the retardance of a 
conventional Cassegrain. 

In general, a Cassegrain telescope displays both these 
polarization aberrations simultaneously, linear polarization 
defocus and linear retardance defocus, with the retardance being 
the larger term. This polarization aberration induces 
polarization coupling which reduces the accuracy of subsequent 
polarimetric measurements. The polarization coupling for a 
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SAMEX-type Cassegrain design but with standard aluminum coatings 
is 3 x 1CT 3 , almost two orders of magnitude greater than the 
radiometric accuracy of the optics. Similar amounts of 
polarization accuracy would be associated with relay optics 
utilizing standard antireflection coatings. 

3.3 The SAMEX Design 

t 

It is clear from the discussion of the previous section 
that standard optical designs will not suffice for the SAMEX 
foreoptics because they produce unacceptable levels of 
instrumental polarization. The design of the SAMEX foreoptics 
resulted from the development of a new method for the analysis of 
instrumental polarization based on the theory of "polarization 
aberrations." This theory allows the description of the 
variations of amplitude, phase, and polarization of an optical 
wavefront across the exit pupil of an optical system. Because 
the theory naturally incorporates the polarization properties of 
the thin film coatings on the individual surfaces of the optical 
system, this method integrates the coating design with the lens 
design. In ordinary optical design work, these two phases of the 
design are normally decoupled and pursued separately. 

This unified optical and coating design was performed 
for the SAMEX foreoptics to insure that the SAMEX magnetograph 
will accurately measure the polarization state of incident 
sunlight. To improve the SAMEX polarization accuracy over 
standard designs, various telescope and lens coatings were 
investigated. Most standard reflection-enhancing coatings for 
telescope mirrors were found to be substantially worse than the 
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Figure 28. Linear retardance effects associated with a 
Cassegrain telescope. The large effect of linear retardance from 
a Cassegrain mirror is shown for (a) a beam of uniform left 
circular polarization and for (b) a beam of uniform linear 
(vertical) polarization. The dominant effect of linear 
retardance is a coupling of linear polarization into circular and 
vice versa. 




bare aluminum. Two coatings were designed which had 
significantly improved polarization performance for a 
telescope. These coatings are described in the section on 
coatings (section 3.6). Figure 29 shows the linear polarization 
(a) and linear retardance (b) associated with a Cassegrain 
telescope when these coatings are utilized and when ordinary 
aluminum coatings are used, as a function of the height of the 
ray entering the telescope, from top to bottom. The amount of 
polarization or retardance for the three coating choices is 
plotted along the y-axis for (1) two aluminum coated mirrors, (2) 
two 0201 coated mirrors, and (3) one Q201 coated mirror and one 
eight-layer enhanced-ref lection coating. The third design (3) 
balances polarization and retardance effects of opposite signs 
between the two mirrors to achieve a polarization performance 
superior to either mirror separately. 

Similar design strategies have been used with the relay 
lenses to achieve significant improvements over conventional 
coatings and to balance the remaining polarization effects. As a 
result, the overall polarization performance for the foreoptics 
has been improved by a factor of ten thousand relative to 
conventional designs. Thus the actual attainment of the 
polarimetric sensitivity determined by the SAMEX requirements 
will presumably be limited only by whatever scattering effects 
are present, and not by the instrumental polarization of the 
optics. 

The method used to achieve this dramatic result - 
polarization aberration theory - is outlined briefly in sections 
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Figure 29. Instrumental polarization effects for rays through a 
Cassegrain telescope. Linear polarization (a) and retardance (b) 
as a function of pupil coordinate in a Cassegrain telescope are 
plotted for three different sets of reflective coatings on the 
primary and secondary mirrors. The three sets of coatings are: 
(1) two aluminum (dot-dash curve), (2) two Q201 coatings (dashed 
curve), and (3) one 0201 and one aluminum (solid curve). The 
0201 interference coating is a 201 coating of bi-layers with high 
(n * 1.52) and low (n: .*.1*45) refractive indices. 


3.4 and 3.5 and expanded upon in Appendix C. A full development 
of the method is given in Chipman (1987). The basic results from 
the method are the second-order polarization aberration 
coefficients which provide a quantitative measure of the 
polarization accuracy of any optical system. Specifically, these 
coefficients determine the parameter A , the polarization 
accuracy, defined as the maximum fraction of light which can be 

i 

coupled into an orthogonal polarization state. It is given in 
terms of the second order polarization aberration coefficients 
P(l,0,2,2), P(l, 1,1,1), and P(l,2,0,0): 

A p < P(l,2,0,0) 2 + P(l,l,l,l) 2 + | P(l,0,2,2) 2 . 

For the SAMEX magnetograph foreoptics and coating design given 
herein, the value of 

Ap < 1.4 x 10 -7 

was obtained. For standard coatings of aluminum, the 
polarization accuracy is 

A = 2.7 x 10" J 
P 

for a Cassegrain telescope alone (no relay lenses). (The second 
order aberration coefficients used in these calculations are 
given in Table 19.) This result for Ap means we have achieved 
our design goal: the polarization state of the light from the 
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Sun can be determined without introducing polarizing effects from 
optical elements in the system. We have been able to effectively 
eliminate the problem of induced instrumental polarization in the 
SAMEX foreoptics. 

3.4 Polarization Aberration Theory 


The SAMEX foreoptics are intended to transmit all 
polarization states equally. But all optical interfaces display 
some polarization when used at non-normal incidence. Thus 
polarization is present in all systems at some level. If the 
system is intended to be nonpolarizing, the instrumental 
polarization is often termed "residual polarization" to signify 
its generally undesirable character. Residual polarization might 
be compared to wavefront aberration because both interfere with 
the measurement of optical fields and reduce the image forming 
potential of the optical system. 

The principal cause of instrumental polarization in the 
optical systems of present solar magnetographs is the 
polarization due to non-normal incidence at the optical 
interfaces and coatings. Since each ray takes a different path 
through the system with its own angles of incidence and planes of 
incidence, each ray in general experiences a different change in 
its state of polarization. This residual polarization varies 
with wavelength, object coordinates and pupil coordinates. 
"Polarization aberrations" will be defined as variations of the 
amplitude, phase and polarization of an optical wavefront across 
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the exit pupil of an optical system and the dependence of these 
variations on wavelength and object coordinate. The polarization 
aberrations are extensions of the wavefront aberrations and 
encompass both amplitude and polarization variations, thus 
providing a more complete characterization of the electromagnetic 
fields transmitted by an optical system. 

Vacuum-deposited thin filings are used on most optical 
surfaces to control the amount of transmission and reflection. 
These thin films are usually less than the wavelength of light in 
thickness. Being so very thin, the effect of the films on ray 
paths are accurately modeled by treating the films as having 
parallel surfaces which contour the substrates on which they are 
deposited. Due to the closely spaced parallel surfaces, thin 
films have negligible influence on the ray paths through the 
system and are generally ignored when simulating a system by ray 
tracing. These coatings principally affect the amplitude and 
polarization of the ray and have much less effect on the optical 
path difference. This division, with the optical surfaces 
governing the ray paths and the thin film coatings governing the 
amplitude and transmission, allows the optical system design 
problem to neatly decouple into two separate problems, lens 
design and coating design. The wavefront performance and image 
quality of the system is calculated by a lens designer using a 
ray tracing optical design program. The amplitude and 
polarization calculations at individual surfaces are performed 
using a thin film design program. 
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This decoupling of optical design and coating design has 
usually worked well. The coatings designed to optimize the 
transmittance or reflectance at an interface have usually reduced 
the amplitude and polarization variations and thus reduced the 
polarization aberrations at the interface as well. For example, 
a quarter-wave magnesium fluoride antireflection coating on glass 
typically reduces reflection losses. at the design wavelength by a 
factor of four, and reduces the instrumental polarization by a 
comparable factor. This fortuitous circumstance has allowed lens 
and coating design to remain relatively uncoupled. Thus 
instrumental polarization was usually ignored as a higher order 
effect. But it is not sufficient to design thin film coatings in 
isolation from the lens design for the SAMEX magnetograph - the 
demands on amplitude and polarization performance are too great. 

For the designs of the SAMEX system special methods 
have been developed to calculate the instrumental polarization of 
the SAMEX foreoptics. These methods are described in detail in 
Appendix C. Calculating the instrumental polarization requires 
performing thin film calculations during the ray tracing 
process. This idea is not new, but its implementation is complex 
enough to have delayed this obvious integration of these two 
branches of optical design until specifications required it. 

In this new methodology, a Jones polarization matrix is 
calculated for arbitrary optical paths through the optical system 
and includes the effects of all the specified optical coatings on 
the curved optical interfaces. The specific technique used for 
this calculation is the method of "polarization aberrations." 
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This is a method analogous to the aberrations of geometrical 
optics (spherical aberration, coma, astigmatism, etc.) except 
that it encompasses the polarization effects of coatings as well 
as the wavefront aberrations. 

The analysis of the SAMEX optics proceeds in several 
stages and is summarized in Table 16. First, the optical system 
is designed using the CODE V optical design program to optimize 
the optical design for a spatial resolution with a minimum 
aperture. During this phase the angles of incidence are kept as 
small as possible to reduce polarization effects from coatings. 
Second, special thin-film, reflection-enhancing coatings are 
designed for the telescope mirrors and special antireflection 
coatings are designed for the lenses because conventional .coating 
designs show significant polarization and retardance 
contributions near normal incidence. These coatings are special 
designs which reduce the polarization effects at small angles of 
incidence over the wavelength band of the magnetograph. Next, a 
Taylor series is calculated to represent the coating 
performance. Finally, the effect of these coatings in the 
magnetograph optical design is calculated using the polarization 
aberration method; this last step produces the parameter the 
polarization accuracy. 
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Table 16. Steps in Instrumental Polarization Simulation 


1. Design optical system using CODE V. 

2. Design low polarization, high reflectivity telescope coatings. 

3. Design low polarization, antiref lective lens coatings. 

4. Determine Taylor series representation of coating 
performance . 

5. Calculate polarization aberration .coefficients of optical 
system . 

6. Calculate A . 

P 

7. Iterate 2, 3, 4, 5- 6 until satisfactory performance is 

achieved (A < 10"°). 

P 


In the following section, the method for deriving the 
polarization aberration coefficients and the polarization 
accuracy A^ is outlined in more detail. 

3.5 Polarization Aberration Coefficients and Polarization 
Accuracy 

The polarization states of the electromagnetic field 
are described by the complex two-component Jones vector, 3. The 
polarization states and hence the Jones vector are transformed 
when a ray passes through an optical interface which is described 
by the Jones matrix transformation JJ. For the SAMEX foreoptics 
we have homogeneous, weakly polarizing optical elements (by 
design) for which the transmission coefficients perpendicular, 
t s , and parallel, tp, are given in terms of a Taylor series 
expansion of the angles of the chief and marginal rays, i c and 
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i f respectively. Each individual optical ray path for the 
rotationally symmetric system can be defined in terms of the 
object height, H, and the pupil coordinants, p and as defined 
in Figure 30. The cascade effect of an optical train is given by 
the product of the individual Jones transformation matrices for 
each element interface. Then the overall polarization of the 
foreoptics is described by the complex Jones transformation 
matrix which is expanded in terms of the ray coordinates 
(p, <j>, and H) and the basis matrix set a(k): 

3 

JJ = E E E E P(k,u,v,w) H u p v cos w U) a(k ) 

k=0 u v w 

where P(k,u,v,w) are the expansion coefficients. The 2 by 2 
matrices, <x(k), which define the k th polarization state, are- the 
identity matrix and the Pauli spin matrices. This expansion is 
in the same form as the standard wavefront aberration expansion. 
These polarization expansion coefficients for an expansion to 
second order (specified by s = 2) in the angles of incidence are 
a function of the total optical transmission, t, the normalized 
secondard order basis set coefficients, d(k,s,q), for each 
individual optical element, q, and the angle of incidence of the 
chief and marginal rays. Hence these coefficients are dependent 
on the characteristics of the optical coatings and the ray 
tracing results. These SAMEX specifications are given in Table 
19a-c . 

The polarization aberration expansion for the radially 
symmetric system of interfaces with isotropic coatings and 
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Figure 30. Paraxial coordinate system. The paraxial system is a 
normalized right-handed coordinate system. The z axis is the 
optical axis of a rotationally symmetric optical system; light 
initially travels in the direction of increasing z. Rays through 
an optical system are characterized by ray coordinates at the 
object and entrance pupil. H is the normalized object 
coordinate/ p is the normalized pupil radius/ and ♦ is the polar 
angle in the pupil measured counterclockwise from the y axis. 

The normalized Cartesian coordinates in the pupil are x and y. 

The chief and marginal rays are also shown. 


nonpolarizing transparent media is given to second order in the 
angles of incidence by the polarization aberration expansion 

JJ(H,p,*) = T[a Q + P(l,2,0 / 0)H 2 a 1 + 

P( 1 / 1 , 1 , 1 ) Hp ( c^cos* - a 2 sin<fr) + 

P( 1 , 0 , 2 ,2 ) p 2 ( <7^0032$ - a 2 sin2$)]. 

The amplitude a^j]^ of the polarization aberration coefficient 
P(i/j#k,l) describes the magnitude of the polarization aberration 
effects while the phase 5^^ of the coefficient describes the 
magnitude of the retardances. Then, given a specific incident 
polarization state, the polarization state of the exit beam can 
be calculated from this matrix. Therefore the amount of 
polarization in the orthogonal state can be determined. This 
orthogonal polarization in the exit beam then determines the 
accuracy of the polarization measurement. 

As an example of the coupling of the optics to an 
orthogonal state of polarization, i . e . , polarization crosstalk, 
consider the following example. The on-axis linear polarization 
and linear retardance of the SAMEX foroptics, i.e., the term 
linear defocus, is described by the fourth term in the 
polarization aberration expansion JJ. The instrumental 
polarization function J^(H,p,^) for linear defocus is then 

( H , p , $ ) = t [ ct q + p d,0, 2,2)p 2 (a^cos2$ - cj 2 sin2$)] 

r 2 t 

= t [ <Jq + ( a 1022 + l5 io22 p ( ff 1 cos2$ - a 2 sin2<fr)J, 
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where x is the amplitude of transmittance of the system down the 
optical axis. The transmittance (x) describes the polarization 
independent reflection and absorption losses associated with the 
ray down the optical axis. P(l,0,2,2) describes the linear 
polarization (aio22^ an( ^ linear retardance ( 5 iq 22^ associated 
with the marginal ray. 

At any given point in the pupil, the eigenenpolarizat ions 

t 

A 

are linearly polarized light, oriented radially, J r , and 

A > 

tangentially J fc . This concept is important to the definition of 
polarization accuracy which is to be defined. We shall now 
calculate for circular and linear incident polarize light to 
illustrate of the effect of crosstalk. 

Maximum coupling occurs when the incident light is 
circularly polarized, since circularly polarized light can always 

A A 

be decomposed into equal components of J r and everywhere in 
the pupil. The coupling is zero in the center of the pupil 
(where the polarization and retardance vanish) and increases to a 
maximum coupling at p 8 1 of 

I c,max( H ' 1 ' < *’ ) = |P(lrO,2,2)| * a 1022 + 6 1022 

at the edge of the pupil. The net fraction of incident 

circularly polarized light coupled into the orthogonal circularly 

polarized state, 1^, is given by the integral over the pupil, 

0 

2 7T X 

ij - ■! / * d* pdp|p 2 P(l,0,2,2)| 2 , 

_ 1 P( 1 , 0,2,2) | 2 
- 3 . 
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If we have a Stokes vector of +V then the amount of polarization 
in the orthogonal state -V can be calculated. If -V is the 
maximum amount of crosstalk, we identify this as the polarization 
accuracy. 

For incident linearly or elliptically polarized light, the 
fraction of coupled intensity is less because the light is not 
composed of equal fractions of eigenstates. The coupling is 
minimum for incident linearly polarized light, which will be in 
one of the eigenpolar izations along one axis in the pupil and in 
the orthogonal eigenpolarization along the orthogonal axis. Here 
the fraction of coupled energy will be calculated assuming an 
incident polarization state of horizontal linearly polarized 

A 

light H for calculational simplicity; the same fraction is 
coupled for any linearly polarized incident state. The 
orthogonal state of vertical linearly polarized light is 

A 

designated as V. The polarization state transmitted by an 

A 

optical system described by linear polarization defocus for H is 


3(H,P,+) = x [H + P( 1,0,2, 2)p 2 (H cos 2$ - V sin 2<fr)]. 


The fraction of incident H light coupled into V light is equal to 


I 2 = i ^ d* / pdp|«?(H,p,$) • v | 2 


1 2 7T 1 2 

= — / d<j> / pdp|p(l,0,2,2)p sin 2<J> 

* 0 0 
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P(l,0,2,2) 

IT 


4, W ^ 

/ sin 2$d$ 
0 


/ 5 P dp 
0 


P{ 1 , 0 , 2 , 2 ) 
g 


Therefore, we defined the polarization accuracy, A , as the 
maximum fraction of light which can be coupled into orthogonal 
polarization states. The incident polarized state is given by 
the Jones vector, 3. This vector is rotated by the optical 
system and the rotation is given by the Jones matrix, JJ. The 
amount of polarization along the orthogonal state of 
polarization, 3', of the incident polarization state is given by 
the projection of JJ(J) into 3'. i.e. 


JJ(3) *3' . 


The maximum projection is used in defining the polarization 

1 2 i 

accuracy. In the above example I_ > I c , and only I c is 
retained for the linear defocus term, P(l,0,2,2). 

This value is given by the square of the second order Jones 
matrix and is given in terms of the polarization aberration 
coefficients. Then the polarization accuracy is given in terms 
of the second order polarization abberation coefficients 
P ( 1 , 0 , 2 , 2 ) , P(l, 1,1,1), and P(l, 2,0,0): 

A = — / I JJ 0 (l,p,<j)) 1 2 pdpd$ 

P IT ‘ 1 Z I 
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or 


A p < |P(1,2,0,0 ) | 2 + \ | P(l,l # l,l ) | 2 + 3 | P( 1,0, 2, 2) | 2 , 

where the integrations have been carried out for the squared 
terms and estimated for the crossed terms as being less than or 
equal to the direct product terms. For the SAMEX magnetograph 
design, the value for the polarization accuracy, 

A < 1.4 x 10 ~ 7 
P 

is obtained for specially selected optical coatings. These 

second order aberration coefficients are given in Table 19. The 

second order coefficients are sufficient since the next order 

that contributes is the fourth order and their polarization 

2 

effects would be on the order of ( A ) . In the following section 

P 

we discuss the process used to design these special optical 
coatings . 
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Table 19— a. Table of Total Second Order Polarization 
Aberration Coefficients for the 
SAMEX Magnetograph Given the Individual 
Surface Components of Table 19-b 


Vector Quadratic Piston : 

P(l,2,0,0) = * d(l,2,q) i c 2 (q) 

= -7.72x10“^ - T 1.50x10“^ radians^ 

Vector Tilt: 

P(l,l,l,l) = 2 * d(l,2,q> i c (q) i m (q) 

= -1.42x10“^ - T 2.70x10”^ radians^ 

Vector Defocus: 

Pd/0,2,2) = “ d( 1, 2,q) i m 2 (q) 

= -1.91x10 ^ + T 1.80x10“^ radians 2 
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Table 19-b. Surface by Surface Polarization Aberration Contributions 
Given Surfaces as Defined in Table 19-c 


Vector Quadratic 

; Piston 

Tilt 

Vector 

Vector 

Defocus 

d ( 1 / 2 ) i c < 
z Re Im 

•(q) 

TTCTTTT 

Re 

i c (q)im(q) 

Im 

dTTTTj 

Re 

V (q) 

Im 


1 

8 . 3e-10 

8 . 4e-7 

-3 . 9e-8 

-3.9e-5 

1 . 8e-6 

1 . 8e-3 

2 . 

-l.le-8 

-1 . 3e-6 

3. 9e-7 

4 . 6e-5 

-1 . 4e-5 

-1 . 7e-3 

3 

-4 . 7e-7 

-l.le-6 

4. 4e-7 

5.8e-7 

-1 . 3e-7 

-3 . Oe-7 

4 

1 . le-5 

-2 . 6e-5 

-3. le-6 

. -7 . 5e-6 

-9 . le-7 

-2 . 2e-6 

5 

-4 . 9e-6 

-1 . 2e-5 

3. 6e-7 

8 . 5e-7 

-2 . 6e-8 

-6. le-8 

6 

-1.3e-6 

0 

-2. 3e-6 

0 

-3 . 8e-7 

0 

7 

-5 . le-6 

-1 . 2e-5 

-1 . 4e-6 

-3 . 4e-6 

-4. le-7 

-9 . 8e-7 

8 

-2 . 8e-5 

-6 . 6e-5 

-2 . 3e-6 

-5 . 5e-6 

-1 . 9e-7 

-4 . 6e-7 

9 

-1 . 3e-5 

-3. le-5 

-1 . 5e-6 

-3 . 6e-6 

-1 . 7e-7 

-4 . 2e-7 

10 

-l.le-8 

l.le-6 

5. le-8 

-4 . 9e-6 

-2 . 3e-7 

-2 . 2e-5 

11 

-5 . 7e-7 

0 

1 . 2e-6 

0 

-2 . 7e-6 

0 

12 

-l.le-6 

-2 . 7e-6 

1 . 3e-6 

3. le-6 

-1 . 5e-6 

-3 . 6e-6 

Note 

: The surface numbers 

S(i) are 

related to the 

surface 

numbers 


q by S ( i ) 

* q + 2 , 

i.e. q=l 

is the primary 

mirror 

and the 


prefilter 

is not considered. 

The notation 

8 . 3e-10 

means 8.3x 
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Table 19-c 


The Surface Definitions for the Polarization Aberration 
Coefficient Calculations 


Surface 

q 

Optical 

Coating 

Paraxial Anqle 

Second Order ^ 
Polarization 
Re Im 

Magnitude o: 
2nd Order 
Polarizatioi 

1 

refl201 

0.00134 

-0.06250 

2.82e-7* 

2.84e-4 

2.8e-4 

2 

layers 8onAl 

0.00236 

-0.08519 

-1 . 20e-6 

-1 . 39e-4 

1.4e-4 

3 

vl0c5256 

0.03862 

-0.02031 

-1 . 90e-7 

-4.53e-7 

4.9e-7 

4 

n 

0.18664 

0.05388 

It 

i 

It 

It 

5 

w 

0.12576 

-0.00908 

ft 

It 

It 

6 

kf9fk5 

-0.30317 

-0.05091 

-8.91e-8 

0.0 

8.9e-8 

7 

vl0c5256 

-0.12743 

-0.03622 

-1.90e-7 

-4 . 53e-7 

4.9e-7 

3 

II 

-0.29808 

-0.02495 

it 

It 

It 

9 

It 

-0.20326 

-0.2363 

n 

If 

It 

10 

vl0c525 

-0.00703 

0.03169 

-1.38e-7 

1 . 33e-5 

1.3e-5 

11 

kf9fk5 

-0.06243 

0.13588 

-8.91e-8 

0.0 

8.9e-8 

12 

vcl0c5256 

-0.5980 

0.06974 

-1.90e-7 

-4 . 53e-7 

4.9e-7 


*(2.S2e-7 = 2.82xl0~ 7 ) 

**The rate of linear polarization (Re) and retardance (Im) per angle of incidence 
squared (deg -2 ) for the specified optical coating. 
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3.6 Coating Designs 

A design study was performed to find coatings especially 
suitable for use in the SAMEX magnetograph. First, many 
conventional coatings were examined and found to have 
polarization and retardance properties that make them 
unacceptable for use in the SAMEX magnetograph. Alternative 
coating designs were devised with improved polarization 
performance and the best of these were specified for the 
instrumental polarization simulation. Coatings were designed, 
optimized and analyzed using two thin film coating programs: FTG 
Software's Filmstar and Prof. A. Macleod's (University of 
Arizona) FILMS. 

The requirements for the SAMEX magnetograph coatings ar.e the 
following : 

1) They must have extremely low polarization properties 
near normal incidence so that they introduce a minimum of 
polarization or retardation into the optical system. 

Specifically, the differences between the s and p amplitude 
coefficients for reflection or transmission for less than 10° 
angle of incidence should be 0.01% and the phase difference 
should be less than 0.05°. 

2) The system's transmission must not be compromised by the 
coatings: reflective coatings should reflect > 99.9% and 
transmission coatings transmit > 99.7%. 

3) They must be manufacturable and not be excessively 
costly . 
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4) The amount of scattered light should be small. 

5) They should have reasonable tolerances for fabrication. 


Several common reflecting coatings for telescopes and 
antireflecting coatings for lenses were analyzed and found to be 
inadequate for the SAMEX magnetograph. For example, telescopes 
are frequently coated with bare aluminum coatings (Al) or 
quarter-wave stacks on top of aluminum. These coatings have very 
little linear polarization near normal incidence but substantial 
amounts of linear retardance. The common antireflection coatings 
for lenses (quarter-wave MgF 2 and broadband antiref lective 
coatings) have very little retardance near normal incidence but 
substantial linear polarization. 

Three coatings (see Table 20) have been selected from our 
coating design study for inclusion in the SAMEX magnetograph 
design. The polarization effect of these coatings has been used 
in the instrumental polarization simulation. These coatings 
demonstrate that coatings with the necessary low polarization 
performance can be designed. The present coatings are somewhat 
sensitive to fabrication errors and are very wavelength 
sensitive. The designs would benefit from further work to 
understand why they work, knowledge which should lead to fully 
optimized designs with improved manufacturing tolerances. 

Table 20 contains the description of the three coatings 
specified. In addition, the optical system contains two cemented 
interfaces, between kf9 glass and fk5 glass, specified kf9fk5. 


The polarization of this interface, although negligible, was 
included in the instrumental polarization calculation. 
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Figure 31. "Optical properties of the interference coating Q201. 
The phase and reflectance properties of the interference coating 
Q201 designed for the SAMEX telescope mirrors are shown as func- 
tions of the angle of incidence (A, on the left) and wavelength 
(B, on the right). (Al) The phase of a reflected beam is shown 
as a function of the angle of incidence. The s and p components 
of an incident wave of wavelength 5250 A are shown as one curve. 
(A2) The difference in phase of the s and p components shown in 
(Al) is plotted as a function of the angle of incidence. (A3) 
The % reflectance of a 5250 A reflected wave is shown as a func- 
tion of the angle of incidence. (Bl) The phase of the reflected 
wave is plotted as a function of wavelength. (B2) The percent 
reflectance of a wave incident at an angle of 2° is shown as a 
function of wavelength. (B3) The percent reflectance of a wave 
incident at an angle of 6° is shown as a function of wavelength. 








The coating Q201 is a 201-layer, quarter-wave enhanced 
reflection coating specified for the primary mirror. The 
philosophy behind this design is that the polarization effects 
arise because of large refractive index differences between 
layers. This coating uses two materials with a small index 
difference, fused silica at n=1.45 and evaporated glass at 
n=1.52. To achieve the desired reflection performance with 


smaller index differences requires more layers, 201 in this 
case. The materials chosen are both amorphous and thus quite 
suitable for coatings with large numbers of layers, since they 
don't display the microcrystalline growth patterns which lead to 
the unacceptable scattering and inhomgeneity associated with 
crystalline materials. Figure 31 shows performance curves for 
Q201 both as a function of angle of incidence and wavelength. 
Figure 31.a-l shows the absolute phase change on reflection for 
the s and p components. The quadratic portion of these curves is 
the defocus introduced by the coating. The parameter of greatest 

iuicLcat i.o Cue i. e laiuan^c , uuc uii.i.3isii^c u>si.»»cdi unc o a Hu y 


phase changes; it is plotted separately in 31.a-2. This 
retardance is below 0.1° over the angles of incidence of 
importance. Figure 31.a-3 shows the s and p intensity 
reflectance as a function of angle of incidence. The difference 
between the s and p reflectance is the linear polarization, here 
less than 0.01%. 

The coating specified for the secondary, L8onAl, is a more 
conventional 8-layer, reflection-enhancing coating specially 
optimized to complement the Q201 coating. This coating has the 
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opposite sign on both its linear polarization and retardance 
relative to Q201. Thus the small residual polarizations of these 
two coatings tend to cancel, resulting in a nearly polarization 
free Cassegrain design. 

The coatings specified for the lens surfaces are a two-layer 
coating of the V coating family of designs. A region of solution 
was found with the remarkable property that the linear 
polarization and linear retardance both changed signs within 50 
nm in wavelength of each other. Thus, by varying the thicknesses 
of the layers, low polarization coatings with several different 
useful properties were obtained. The two scalings used in the 
design were 1.0 for the coating vl0c525 and 1.06 for vcl05256. 
These coatings have by far the lowest polarization effects of any 
lens coatings investigated. Fabricating and testing samples of 
these coatings should be conducted in the next study phase. 


Table 20. Coatings Specified for the Magnetograph. 


Name 

Type 

Number of 
Layers 

Polarization 

Magnitude 

Reflectance or 
Transmission 

Coating 

Design 

V10C525 

Antireflection 

2 

4.9e-7 

99.7% 

( ,25H ,25L)100 .254H 
H=1.52 L=1.45 

Q201 

.82L) 


Reflective 201 

2.8e-4 

99.98% ( .062H 
H=2.15 L=1.38 

L8onAl 
. 25L)2 

Reflective 

9 

1.4e-4 

99.65% 

.38H .29L .26H 

.25L ( .25H 

H=2.38 L=1 .38 
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All the polarization analyses of thin films performed here 
have assumed ideal/ thin-film structures, materials which are 
uniform, isotropic, homogeneous and free from scattering. This 
is a sufficient assumption for the design and analysis of the 
thin films. Actual thin films possess extremely complex 
microstructure which contributes to scattering, inhomogeneity and 
anisotropy. The typical thin film microstructure is an array of 
columns growing up out of the substrate. This causes nonideal 
coating properties which will impair the performance of the 
system at some level. 

The performance of the SAMEX magnetograph will be enhanced 
by using the finest coatings available, which will probably be 
deposited using ion-assisted thin film deposition. By bombarding 
the growing thin film with energetic ions, usually helium or 
argon, the growing thin film structure is pressed down and 
becomes denser (Martin, 1986). This disrupts the growth of the 
large columnar structures which are responsible for much of the 
scattering and anisotropy. The resulting films have lower 
scatter and greater uniformity. The ion-assisted thin film 
deposition technique is still new and not yet widespread in its 
use . 

As part of the preliminary efforts leading to the actual 
fabrication of the SAMEX magnetograph, test coatings need to be 
prepared and coating vendors qualified. Several coating 
manufacturers should be contracted to produce small (2x2 inch) 
samples of the specified coatings. These coatings should be 
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tested for spectral transmission, polarization performance and 
scatter (bidirectional reflection measurements) as part of the 
final design and vendor qualification. The Marshall Space Flight 
Center and the University of Alabama in Huntsville Center for 
Applied Optics have the expertise for such coating 
characterization. These data would then be used to calculate the 
impact of the stray light coming from the coating scatter on the 
magnetic field accuracy of the magnetograph. 

3.7 Further Sources of Instrumental Polarization 

Two effects that could degrade the polarization performance 
of the SAMEX optics are coating anistropy and scattered light. 
Coating anisotropy is the variation of refractive index with 
direction in the coating. Coating anisotropy can be measured 
ell ipsometrically . It frequently occurs where coatings have been 
deposited at non-normal incidence. Most coatings have a columnar 
microstructure which usually grows out of the substrate toward 
the source. For a nonnormal deposition angle, the resulting 
coatings have builtin birefringence. For the magnetograph 
optical design, this additional coating-induced retardance is a 
problem which must be held to acceptable levels in the final 
coatings. A level of anisotropy below a root sum squared 
birefringence of 0.0001 radians rms per coating is highly 
desirable . 

Scattering depolarizes light: the scattered light is random 
and carries less information about its original polarization 
state. Coatings, because of their detailed microstructure, may 
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display substantial scattering. Considerable effort is being 
devoted on many approaches to producing optics with reduced 
scatter/ including the ion-assisted deposition coatings already 
mentioned. Scattering in the coatings of the foreoptics will 
probably be the limiting factor in the accuracy of the 
magnetograph/ now that the instrumental polarization of the 
Cassegrain telescope has been reduced by orders of magnitude. 

3.8 Optical Tolerances 

An optical tolerance study of the magnetograph is outside 
the scope of this study. However we would like to point out the 
areas of special study that are needed for a flight instrument. 

The individual optical components need to be qualified for 
spaceflight. This would include operation under a vacuum for an 
extended period of time. The effects of particle radiation over 
the lifetime of the spacecraft must be resolved. The mechanical 
robustness of the optical system is to be defined. This would 
include the effect of thermal drift on the optical alignment. 

The general effects of mechanical fitting and sensitivity due to 
fabrication errors need to be addressed, including surface 
figure, tilts, decentrat ions , and optical coating variations. 

The overall alignment and calibration procedures are to be 
specified before the optical design is completed (Yoder, 1986). 


31 


APPENDIX C 


POLARIZATION ABERRATION THEORY 

Introduction 

This appendix describe the method which have been developed 
to calculate the instrumental polarization of the SAMEX optics. 

A full development of the method is contained in "Polarization 
Aberrations" by Chipman (1987). The calculation requires 
performing thin film calculations during the optical design 
process to determine the two parts (eigenvectors) of an optical 
beam, as a function of the object and pupil coodinates. 

The mathematical method of treating instrumental 
polarization will be first discussed in terms of the Jones and C- 
vectors and then the polarization aberration expansion will be 
derived . 

The Jones Matrix and C Vector for the Characterization of 
Polarization 

The most efficient mathematical method for treating the 
SAMEX instrumental polarization is the Jones calculus. The 
Mueller calculus is a more difficult representation which 
includes optical depolarization (scatters) properties. Such weak 
depolarization effects are more readily handled by experimential 
measurements. The Jones calculus (Jones 1941a-c, 1942, 1947, Azzam 
and Bashara 1977, Theocaris and Gdoutos) is a mathematical 
formalism for treating problems involving the description of 
polarized light and polarizers which uses the Jones vector for 


the description of polarized light and the Jones matrix to 
characterize the polarizing properties of an optical element. 

The details of the following discussion are given by Chipman 
(1987). The elements of the Jones matrix and the C vector 
characterization are outlined first. 

Definition of the Jones Vector in terms of the Electric Field 
Amplitudes 

The Jones vector expression for the a quasi -monochroma tic 
plane wave propagating parallel to the z axis with electric field 
amplitude 

E(t) = E x (t) + E y (t) , 


where, 

E x (t) = n x E 0,x cos [ (kz-wt ) + e x l 


and , 


E y (t) = n y E 0ry (t) cos [(kz-wt ) + e y ] , 

where n x and n y are direction unit vectors in the x and y 
direction with the light propagation in the z direction. There 
are 4 parameters Eq, x , E o,y' e x' anc5 e y defining the wave beside 
the wavelength. 

The time dependent Jones vector is defined in terms of the 


electric field amplitudes as, 


The components of J(t) are the instantaneous components of E(t). 

The normalized Jones vector J is a time independent 
normalized vector where all the vector components of J(t) have 
been divided by the incident electric field amplitude, 

J = J(t)/E 0 (t) . 

The normalized Jones vector is referred to as "the Jones vector* 
unless otherwise stated. Knowledge of J and E Q provides all the 
information necessary to reconstruct E(t) to within a constant 
phase factor. 

Table C-l lists the Jones vectors for the most common 
polarization states: horizontal linear , * vertical linear, +45 
degrees linear, -45 degrees linear, right circular and left 
circular polarized light . These vectors can be multiplied by an 
arbitrary phase factor without changing the polarization ellipse 
of the light; it only changes the absolute phase. 

Having established the vector which defines the polarization 
state we now consider the matrices which represent the 
polarization effect of the optical elements and allow 
polarization calculations to be performed. 



M 


Table C-l. The Basic Jones Vector Representation for Linear 
and Circular Polarized Light 



Linear Polarized 

Light 


Horizontal Vertical +45 Degrees 

-45 Degrees 


1 

1 0 1 s 

1 1 

si 1 1 


0 

1 1 ! 

1 | 

l-l 1 


Circular Polarized Liaht 



Right Circular 

Left Circular 


s| 1| 

s 

1 1 


l-i i 


li 1 


J2 

where s = — =■ 
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In his original paper, Jones (1941a) shows that the 
relationship between the Jones vector incident on a polarizer, J, 
and that transmitted or reflected by a polarization element, J*, 
can always be related by a matrix, the Jones matrix, JJ. Only 
certain transformations of the field components are allowed, 
those describable by a matrix. Thus the fundamental relationship 
between the vector components of the electromagnetic fields 
before and after a polarizing element is. 


J' = JJ J 


The Jones matrix, JJ, is a two by two matrix with complex 
elements , 




JJ 


j(l,l) 

j (2,1) 


j ( 1 / 2 ) 
j ( 2, 2 ) 




where j(k,l)*a(k,l)+ib(k,l). 

Thus the Jones matrix has eight degrees of freedom. Thus 
there are eight different forms of polarization behavior. These 
eight forms are listed in Table 2E. Every Jones matrix 
corresponds to a physically realizable polarizer. 

Sequences of Polarizers 

The Jones matrix associated with an optical ray path through 
a sequence of polarization elements is just the matrix product of 
the Jones matrices for the individual polarizers. If an optical 
ray traverses a series of elements, 1, 2, ... Q, and the Jones 
matrices appropriate to that ray for each element are, JJ(1), 
JJ(2), ... JJ(Q), then the Jones matrix describing the 
polarization properties of the system along this ray path is 
given by the matrix product, 

JJ = JJ(Q) ... JJ ( 2 ) JJ(1) 

Since the Jones matrix of an optical element is dependent 
upon the wavelength, angle of incidence, orientation, and path 
through the element, each ray in each wavelength will usually 
have a different Jones matrix. Only if a collimated 
monochromatic beam through a series of planar optical interfaces 
can be assumed, then a single Jones matrix can be written for the 
entire cross section of the beam. This is the case in the 
polarimeter section of the magnetograph but not in the 
foreoptics. 


Coordinate System 

The Jones matrix is defined relative to an arbitrary x and y 
coordinate system. Since these coordinates have been defined for 
the Jones vector/ the coordinate system of the Jones matrix is 
defined in terms of the Jones vector coordinates. 

It is often desirable to align the Jones vectors coordinates 
with the s and p planes of an optical interface. Only for plane 
surfaces does the orientation of the s and p planes remain fixed 
across the surface. For nonplanar surfaces, it is necessary to 
maintain two sets of coordinates, the global x and y coordinates 
with respect to which the Jones matrix is defined, and a local s 
and p coordinate for each individual point on the interface. The 
local s and p coordinate system will have its x' and y' axes 
aligned with the local s and p planes of the surface for the 
evaluation of Jones matrices at given surface coordinates. Then, 
these local Jones matrices will be rotated to bring all the 
matrices into the global coordinate system. 

Pauli Spin Matrix Basis and the C Vector 

The Pauli spin matrices form a most useful basis for 
interpreting the Jones matrix "space" and define a basis set for 
the JJ matrix. The identity matrix, cr ( 0 ) and the Pauli spin 
matrices, cr ( 1 ) , a(2), and a(3), are defined by: 


cr ( 0 ) 


o(l) 

0(2) 

o(3) 

1 0 


1 0 


0 1 


0 -i 

0 1 


0 -1 


1 0 


i 0 



— 


— 


— — 


An arbitrary Jones matrix will be expressed as, 

3 

J J = z c ( k ) <J ( k ) . 
k=0 

The c's are formed into a four element complex vector, the 
"C vector". The C vector expression 

C = [ c( 0 ) , c(l), c ( 2 ) , c{ 3 ) ] , 


an equivalent representation of the Jones matrix. 


JJ = c ( 0 ) o( 0) + c(l) o(2) + c(2) <j(2) + c(3) a( 3) 


c ( 0 ) + c(l) 


c ( 2) - ic ( 3 ) 


c ( 2 ) + ic ( 3 ) c(0) - c(l) 


When necessary, p and $ refer to the amplitude and phase portions 
of the C vector elements. 


C = [p ( 0 ) exp(7$(0) ) , p ( 1 )exp( T${ 1 ) ) , p(2)exp(7*(2) ) , p ( 3 ) ) exp ( Tp ( 3 ) ) ] 

The elements of C are related to the Jones matrix elements by the 
equations : 


s 


2T) . 

The elements of the Jones matrix are related to the elements of C 
by the equations: 

j(l,l) = c ( 0 ) + c ( 1 ) , j(l,2) = c ( 2 ) - 7 c(l) , 

j ( 2, 1 ) = c( 2) + T c ( 3 ) , j ( 2 , 2 ) = c(0) - c(l) . 

The C vector, like the Jones matrix, has eight degrees of 
freedom. Table C-2 contains a description of the meaning of the 
real and imaginary parts of the c vector elements. 


c ( 0 ) = (j(l,l) + j ( 2 , 2 ) )/2 , e(l) = ( j ( 1 , 1 )- j ( 2, 2 ) )/2 

c(l) = ( j ( 1 , 2 ) + j (2, l))/2 , c ( 3 ) = ( j ( 1 » 2 )- j<2, !))/(- 


TABLE C-2. Interpretation of the C Vector Elements 


Matrix 

Coefficient 

Meaning 


o(0) 

p(0) 

Amplitude 

Absorption 

o(0) 

♦ (0) 

Phase 

Phase 

0(1) 

P(l) 

Amplitude 

Linear Polarization along Axes 

0(1) 

4* ( 1 ) 

Phase 

Linear Retardance along Axes 

o(2) 

P ( 2 ) 

Amplitude 

Linear Polarization, 45 deg 

o(2) 

M2) 

Phase 

Linear Retardance, 45 deg 

o(3) 

P ( 3 ) 

Amplitude 

Circular Polarization 

o(3) 

M3) 

Phase 

Circular Retardance 


The Jones Matrix and C Vectors for Specific Polarizers 

Tables of Jones matrices for various polarizers are found in 
Azzam and Bashara (1977, Section 2.2.3), Hecht and Zajac (1974, 
Table 8.6), Shurcliff (1962, Appendix 2), and Theocaris and 
Gdoutos (1979, Table 4.1). A table listing of the Jones matrices 
and C vectors for the most common polarizers and retarders is 
given in Chipman (1987, Table 6). 

The Meaning of the Coefficients of the C Vector 

The primary reason for the introduction of the C vector is 
to simplify the representation of polarizers. Each of the 
elements of the C vector represents a specific type of polarizer 
behavior. 

The real parts of the C vector all correspond to amplitude 
effects, absorption and dichroism. The phase portion of the C 
vector represent phase effects, propagation and birefringence. 

The first element, c ( 0 ) = p ( 0 )exp ( i <fc( 0 ) ) , is the coefficient of the 
identity matrix. Thus it must represent effects that are 
polarization state independent; these are amplitude and phase. 

The last element, c( 3 ) = p( 3 )exp( i $( 3 ) ) , multiplies the spin matrix 
o(3) which is rotation invariant. Thus the c(3) term represents 
the circular polarization effects; p(3) describes circular 
polarization or circular dichroism and <fr(3) describes circular 
retardance or circular birefringence. The remaining two 
elements, c(l) and c(2), represent linear polarization and linear 
retardance. Linear terms require two degrees of freedom: 
magnitude and orientation. Thus, p(l) and p(2) characterize 


linear polarization or linear dichroism, p(l) in the 0 degrees 
and 90 degrees directions, p(2) in the + and -45 degrees 
directions. Likewise, <fr(l) and $(2) characterize linear 
retardance or linear birefringence. 

Rotated Polarizers 

If a polarizer with Jones matrix JJ is rotated through an 
angle 9 (positive if counterclockwise), the Jones matrix becomes 

JJ* (9) = R( 9 ) JJ R(-9). 

The R(9)'s are the Jones rotation matrices: 

R( 9 ) = 

The Jones rotation matrices obey the relations, 

R(a) R(b) = R(b) R(a) = R(a+b) 

and, 

R( a ) R( -a ) = a ( 0 ) . 

The identity matrix is invariant under rotation; 

R ( 9 ) <J(0) R(-9) = a( 0 ) . 

Under rotation, a(l) and a(2) couple into each other; 


cos ( 9 ) -sin(9) 

sin(9) cos(9) 



R( 9 ) a ( 1 ) R(-9) = cr( 1 )cos ( 29 ) + cr(2)sin(29) 

R( 9 ) a ( 2 ) R( -9) = -a( 2)sin( 29) + <j( 2) cos ( 29) . 


o(3) is invariant under rotation; R(9) a(3) R(-9) = o(3). 

Having established the polarization calculus which describes the 

polarization optics, we now apply these matrices to the SAMEX 

» 

optical system. 


Instrumental Polarization 

Two types of polarization calculations can be performed for 
the SAMEX magnetograph: instrumental polarization and transmitted 
light polarization. The first is the calculation of the 
polarization associated with ray paths through an optical 
system. This determines the instrumental polarization function 
as a function of pupil coordinates for a specified object 
point. The other type of calculation determines the state of 
polarization, such as a Jones vector, transmitted by the system 
along a given ray path for a specified input polarization 
state. By iterating this process, the Jones vector as a function 
of position in the exit pupil is calculated. This Appendix deals 
only with the instrumental polarization calculation. Once the 
instrumental polarization function for the system is known, the 
transmitted Jones vectors are readily determined for all input 
polarization states. 


Polarizers are optical elements which divide an optical beam 
into two parts (Jones vector) and transmit those parts with a 
different transmission coefficient and a different phase. The 
two parts of the beam are referred to the eigenvectors or by the 
more descriptive term, "eigenpolarizat ions . " The two 

eigenpolarizat ions are orthogonally polarized and are transmitted 

» 

by the polarizer with no alteration of their polarization states; 
only the intensity and phase changes. 

The word polarizer will be used to refer to both polarizers, 
such as dichroic or prism type, which have a different 
transmittance for the two eigenpolarizations , and retarders, 
which have equal transmittance but a different phase change for 
the polarizations. 'Polarized Light' by Shurcl if f ( 1962 ) is the 
standard reference on the types of polarizers, their definitions, 
parameters and properties. 

Transparent Systems 

The SAMEX foreoptics are highly transparent and weakly 
polarizing and the following calculational method is optimized 
for this case. The ideal Jones matrix for a ray through a 
transparent nonpolarizing system is 

JJ( ideal) = exp (Td) 

where d is the optical path length for the ray in radians. (The 
complex value (-1 )V 2 is denote by 7 ). The Jones matrix operation 
on the Jones vector, which is composed of the two orthogonal 



amplitude components, define the transmitted state of the ray. 

The above ideal Jones matrix is the identity matrix, which 

signifies that the system has no absorption or polarization. 

Since this is the desired form of the Jones matrix for the SAMEX 

foreoptics, the approach developed here obtains the instrumental 

polarization function as a Taylor series of the system Jones 

* 

matrix in the ray coordinates about JJ( ideal). This approach is 
easily modified for systems which are not highly transparent or 
which contain strong polarizers by performing the Taylor series 
about the Jones matrix for the ray down the optical axis. This 
work deals primarily with this simpler version of the problem, 
transparent systems, to streamline the SAMEX calculations. 

S-P Coordinates 

To handle problems involving light at nonnormal incidence at 
curved optical surfaces, it is necessary to maintain two separate 
coordinate systems: x-y coordinates and s-p coordinates. The x-y 
coordinates are the global x, y and z coordinate system used to 
describe the optical system. The optical axis of radially 
symmetric optical systems coincides with the z axis. 

The s-p coordinates are used to perform polarization 
calculations with the SAMEX coatings which are angle of incidence 
dependent. Most frequently, the functional form of the interface 
polarization is known in s-p coordinates. Thus, the rational for 
using s-p coordinates is that, typically, the Jones matrix for a 
ray at an optical interface will be calculated in the s-p 
coordinates. Then it will be rotated into the system x-y 


coordinates. Once all the Jones matrices for the ray at all 
surfaces have been rotated into x-y coordinates, they can be 
cascaded to give the instrumental polarization along that ray 
path in the system x-y coordinates. 

The s-p coordinates are based on the concept of the s and p 
planes. Consider light with unit wave vector k (normalized to 
one) incident at a surface with normal n. The plane of 
incidence, or "p-plane" is the plane which contains k and n. The 
plane perpendicular to the plane of incidence which contains k is 
the "s-plane". Two unit vectors perpendicular to k are defined 
to form an orthonormal basis for this coordinate system. 


Instrumental Polarization 

All optical elements display variation of their Jones 
matrices as the angle of incidence changes. Further, this always 
involves more than just a variation in the intensity and phase of 
the light; it also involves polarization and retardance. A fine 
optical element used in a transparent system will not display 
polarization effects at normal incidence; it may show some 
absorption, reflection loss or phase shift, but not polarization 
or retardance. 

The Jones matrix can be decomposed and can be expressed as 


3 

= I c(k) a( k) , 
k = 0 


1 0 


1 0 


0 1 


0 -i 


- c 0 ( ) +c l( ) + c 2 ( ) + c 3 ) 

U 01 1 0 -1 10 J i 0 


JJ 


where <j(k) are the identity matrix and the Pauli spin matrices 
which describe the specific polarization state. The c(k) defined 
a vector which then described the polarization properties of the 
element. At an angle of incidence i, the C vector will have the 
form 


(c(0,i) ,c(l, i) ,c(2,i) ,c(3,i) ) = 
(p(0,i)e r *<°' i >, pd.ije 7 ** 1 '^, 


p( 2 , i ) e~^ 2 ' 1 ^ / p(3,i)e T * (3 ' i} 


where each component has an amplitude (p) and a phase ( <fr) part. 
The functional dependences of the C vector coefficients are 
calculated from the Fresnel equations and coating equations for 
the interface. 


Weak Polarizers 

A weak polarization element is defined as a polarizer having 
a C vector such that for some range of i: 

y 2 

| c ( 0 , i ) | 2 » ( | c ( 1 , i ) | 2 + | c ( 2 , i ) | 2 + |c(3,i)| 2 ) . 

Such a polarization element transmits light in a polarization 
state similar to the incident state with only weak coupling into 
other polarization states. The polarization behavior is 
dominated by transmission with only traces of polarization or 
retardance. Any polarization present is at the few percent level 
or less, such that any linearly polarized incident beam has a 
transmission coefficient which varies a few percent or less with 


orientation. Similarly, the retardation is degrees or less, far 
less than a quarter wave ( tt/ 2 ) retarder. Near normal incidence, 
metals in reflection (e.g. telescope mirrors) and dielectric 
refracting interfaces (e.g. relay lens) are weak polarizers. In 
addition, near normal incidence, anti-reflection coated lenses 
used in transmission and metals with reflection enhancing 
coatings are typically weak polarizers for wavelengths near the 
thin film design wavelength. 

Amplitude Transmission Relations 

The amplitude transmission equations for an interface are 
the equations which relate the amplitude and phase of the 
electric fields, E, at an interface. The most general amplitude 
transmission equations for a nonscattering linear interface are: 


“ E’(S)‘ 


a ( ss ) 

E ( s ) 

+ a(ps) 

E(p) ' 

E* (p) 

, 

a ( sp ) 

E(S) 

+ a(pp) 

E(p) 

« 


where for this section, the plane of incidence will be aligned 
with the y axis. This equation is equivalent to the Jones matrix 
equation. 


E'(S) ^ 


a(ss') 

a ( sp M 

E(sT 

E ' ( p ) 


a(ps) 

a(pp) Jj 

E(p) 


For interfaces whose eigenpolarizations are linear polarized 
light oriented parallel and perpendicular to the plane of 



\1 


incidence, the transfer of energy across the interface is 
separable into two uncoupled components which can be written in 
the form: 


E'(s) = a(s ) E(s) = p(s) exp ( i $(s } ) E(s) 

E'(p) = a(p) E(p) = p(p) exp( i $(p) ) E(p). 

The amplitude transmission coefficients a(s) and a(p), or in 
polar coordinates, p(s), 4(s), p(p), and $(p), are determined by 
the Fresnel equations for the interface. The type of energy 
transfer equation, where the s and p equations are separable, is 
a "separable amplitude transmission relation." Only polarization 
elements with linearly polarized light as the eigenpolarizat ions 
have the energy transfer equations in the separable form. This 
includes all the elements and coatings used in the SAMEX 
foreoptics . 

The separable amplitude transmission relations correspond to 
a diagonal Jones matrix in s-p coordinates. The Jones matrix and 
C vector for an amplitude transmission interface in s-p 
coordinates are: 


JJ( i ) 


a(s,i ) 0 

0 a (p, i ) j 


and , 


c = 1/2 [ a(s,i)+a(p,i) , a(s, i)-a(p,i) , 0 , 0 ]. 
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Taylor Series Representation of SAMEX Coatings 

In optical aberration theory, expressions for the optical 
path length of ray segments through the the optical system are 
obtained by performing a Taylor series expansion on Snells law, 
the law of reflection and the grating equation, to obtain 
expressions for the optical path length as a power series 
expansion in the ray coordinates. Thus Snells law, 

n sin i = n* sin i' , 
is rewritten for i' as, 

i' = arcsin[ (n/n') sin i ] 
or 

i' = (n/n') i + [(n/n') 3 ) - (n/n’)]i 3 /5 + 0 ( i 5 ) . 

The polarization aberrations are generated in an analogous 
fashion. To obtain the variation of the Jones matrix in the exit 
pupil of a system, the appropriate coating equations are required 
in Taylor series form. For radially symmetric optical systems, 
expansions in the angle of incidence about normal incidence are 
used . 

An isotropic interface appears unchanged as it is rotated 
about the surface normal. Ideally, and to a first approximation, 
the SAMEX coatings are isotropic. So, for any isotropic 


interfaces the Fresnel equations are even functions since the 
surface, does not distinguish between angles of incidence of ,+i 
and -i . 

An even function contains only even terms in its Taylor 

series expansion about the origin. Thus, the Taylor series 

representations of the coating equations has the form, where f(i) 

is a reflection or transmission coefficient, of 

» 

f ( i ) * f 0 + f 2 i 2 + £4 i 4 + 


where 


f n = 1/nl 


d n f ( i ) 
di n 


i=0 ' 


For weakly polarizing interfaces described by amplitude 
transmittance relations, the Taylor series forms of the Jones 
matrix and C vector are calculated as follows. First, the Taylor 
series expansion is determined for the amplitude transmission- 
relations : 


a(s,i) = ag(s) + a 2 (s) i 2 + a 4 (s) i 4 + ..., 
a(p,i) = a 0 (p) + a 2 (p) i^ + ^4(p) i 4 + •••> 


Then, the Taylor series expansion about i =0 in s-p coordinates 
for the Jones matrix is. 


& 


JJ(i) 


a 0 (s) 0 a 2 (s) 

i) = + i z 


0 a 0 (p) 


0 a 2 (p) 


The corresponding C vector expansion in s-p coordinates is 

C = C c(0,0)+c(0,2)i 2 + ... , c(l,0)+c(l,2)i 2 + ... , 0, 0] 

» 

where the nth order c-vector component is given by 

c(0,n) = 1/2 (a n (s) +a n (s)) 

c(l,n) * 1/2 (a n (s) - a n (s) ) . 

A matrix equation to calculate c{0,n) and c(l,n) from thin film 
program results is given by Chipman (1987). For the SAMEX 
coatings characterized by separable amplitude transmission 
relations, the diagonal and circular polarization components, 
c(2,n) and c(3,n), are always zero. The normalized C vector for 
the separable amplitude transmission relations is 


C = t [ l+d( 0 , 2 ) i 2 + ... , d(l,0)+d(l,2)i 2 + ... , 0, 0 ] 


where 

t= c(0,0) and d(k,n) = c(k,n)/c(0,0) 

The Jones matrix and C vector for coordinates other than the 
s-p coordinates are obtained from the polarization rotation 
operation. For example, the s-p coordinates are rotated with 


A\ 

respect to the x-y coordinates by 0, the orientation of the plane 
of incidence. The Jones matrix in x-y coordinates JJ(x,y) is 
related to the Jones matrix in s-p coordinates JJ(sp) by the 
equation 

JJ( x,y ) = R(-9) JJ(s,p) R( 0 ) . 

The Taylor series coefficients for the Fresnel equations 
which govern an uncoated dielectric or metal surface have been 
determined for use in determining the instrumental polarization 
of a conventional Cassegrain telescope. The notation a(s) and 
a(p) will refer to either the reflected or transmitted amplitude 
transmission coefficient, while t(s), t(p), r(s) and r(p) are 
used to refer unambiguously to the transmitted or reflected 
components. The Fresnel amplitude transmission equations are: 

t(s,i) = ( 2 cos i sin i ' )/(sin( i + i ’ ) ) 

= ( 2n cos i)/(n cos i + n'cos i') 

t(p,i) = ( 2 cos i sin i ' )/( ->in( i+i ' ) cos(i-i')) 

= ( 2n cos i)/( n' cos i+ n cos i ' ) 

r(s,i) = (-sin( i-i ' ) )/( sin(i+i') ) 

= ( n cos i - n' cos i')/( n cos i + n' cos i') 

r(p,i) = ( tan( i-i ' ) )/( tan( i+i ' ) ) 

= ( n' cos i - n cos i')/( n' cos i +n cos i') 


The Fresnel equations depend on the ratio of the indices, n and 
n 1 , but not on the values of the refractive indices 
individually. This relative refractive index ratio is defined as 

N = n/n ' 

» 

The Fresnel equations are equally valid for real n, corresponding 
to transparent media, or complex n, corresponding to absorbing 
media or metals. 

The second order Taylor series expansions for the Fresnel 
amplitude coefficients about i=0 are: 


t(s, i ) 

= (2 N)/(N+1) 

+ 

i 2 

( N ( N-l ) ) / ( N+l ) , 

t(p,i) 

= (2 N )/ (N+l ) 

+ 

i 2 

( N 2 ( N-l ) )/( N+l ) , 

r(s,i) 

= ( N— 1 )/ (N+l ) 

- 

i 2 

( N“1 ) / ( N (N+l ) ) , 

r(p,i) 

= ( N-l )/ ( N+l ) 

+ 

i 2 

(N-l )/(N (N+l ) ) . 


The direct method for calculating the Taylor series 
coefficients of a coating series given in the last section are 
impractical for multilayer coatings due to the complexity of 
calculating the partial derivatives of the appropriate amplitude 
transmission equations. The Taylor series coefficients can be 
obtained numerically from the s and p amplitude transmissions 
evaluated at a series of angles of incidence. An algorithm to 
sixth order has been given by Chipman (1987). The algorithm was 
used with the thin film design programs Filmstar and Films to 


obtain the Taylor series expansions of the transmitted and 
reflective amplitudes as a function of the angle of incidence for 
use in the polarization aberration calculations for SAMEX. 

CASCADED WEAK POLARIZERS 

In this section the Jones matrix describing the instrumental 
polarization for light propagating along a ray path through the 
SAMEX foreoptics is derived. Results are also given for the 
instrumental polarization associated with paraxial rays as 
functions of the Taylor series of the C vectors representing the 
optical interfaces. The notation used in this section is 
compiled in Table C-3. 


Table C-3. Notation for Section 


C 


~ f U \ 
\ r\ / 


d ( k ) 
i 

JJ 

k 

d 

q 

0 

v 

r(k) 


a 

t 


p(k) 


C vector 


^ f \ 
U \ /V / 


^ C ■! l' VA i. ^ 1 A L Ak 

tvcaucu 1UUU 


I »• K -i V v** » 

i j_ w x. ui.ai.jf 
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incidence 

Normalized C vector components in s-p coordinates 
Angle of incidence 
Jones matrix 

Pauli spin matrix index: 0, 1,2,3 

Length of a ray segment 

Surface index 

Total number of surfaces 

Orientation of the plane of incidence 

Absorption or polarization coefficient 

Pauli spin matrix 

Normal transmittance 

Phase or retardance coefficient 


Subscript Ordering: k, 1, q. 

For example, d(l,2,3), is the coefficient for: 
the o(l) polarization basis state, that is second 
order in the angle of incidence Taylor series, 
i , for q=3 , the third interface. 


Consider an optical system with Q optical interfaces 
numbered in the order encountered from q=l to Q. No symmetry 
regarding the optical configuration is assumed. Light propagates 
along a specified ray path such as would be calculated by an 
optical ray trace calculation. At each interface some 
polarization is introduced due to differences in the optical 
constants across the interface. In addition, polarization is 
associated with the ray path between interfaces due to optically 
active crystals, dichroism, birefringence, gradient index 
materials or other polarizing mechanisms. But for the 
polarization analysis for SAMEX foreoptics, polarization 
associated with the optical path between interfaces was zero. 
Therefore only interface induced polarization is considered here. 

Homogeneous Optical Systems 

A homogeneous interface has optical properties independent 
of spatial coordinates on the interface. The Jones matrices are 
functions only of the angle of incidence, plane of incidence, and 
optical properties of the interface media, JJ = JJ(i,8,n,n') and 
similarly C = C(i,0,n,n'). The foreoptics and polarimeter 
sections of the SAMEX Magnetograph are homogeneous optical 
systems . 

Likewise, a homogeneous medium has optical properties 
independent of spatial coordinates. An anisotropic crystalline 
medium is homogeneous if it consists of a single crystal. The 
refractive index varies with direction but not with position. 



Radially Symmetric Systems of Lenses, Mirrors and Coatings 

The polarization properties of optical systems comprised of 
lenses, mirrors and "fine" coatings will be developed. A 
radially symmetric optical system has an axis of symmetry, the 
optical axis. It is assumed that the optical elements and 
materials used in transmission are highly transparent and 

i 

nonpolarizing, as is usual in lenses. The polarization 
contribution from the path lengths through highly transparent 
elements is small relative to the polarization arising at the 
interfaces and is neglected. 

The polarization associated with ray paths near the optical 
axis, or in the paraxial regime, will be derived. For this 
paraxial development to be accurate, it is only necessary that 
the angles of incidence are small enough that the polarization 
associated with the interfaces is adequately approximated by a 
second order expansion of the C vector as a function of the angle 
of incidence. For an uncoated lens or mirror, this approximation 
is generally valid for i < 30 degrees. Calculation of the fourth 
and higher order coefficients allows estimation of the accuracy 
of these second order equations. The paraxial region for this 
polarization analysis is typically orders of magnitude larger 
than the paraxial region of geometrical optics (the region where 
the fourth and higher order wavefront aberrations are 
negligible . ) 

Hompgeneous and isotropic interfaces do not display 
polarization at normal incidence. There is only an amplitude and 


Homogeneous and isotropic interfaces do not display 
polarization at normal incidence. There is only an amplitude' and 
phase change which is represented by the complex number, t, the 
normal amplitude transmittance. An isotropic interface such as a 
lens, mirror or coating has a C vector Taylor series in s-p 
coordinates (9 = 0) of the form 

c(i ' 9 ) l 0 * o " t t 1 ' 0 ' 0 ' 0 ! + i 2 t [d(0,2),d(l,2),0,0] + 

+ i 4 t [d(0,4) ,d(l,4) ,0,0] + ... 

For an arbitrary orientation 9 of the plane of incidence, the C 
vector is 

C( i, 8) = t [1,0,0,0] + i 2 t [c(0,2)t,c(l,2)t,c(2,2)t,0r 

+ i^ t [c ( 0 ,4 ) t , c( 1 , 4 ) t , c ( 2 , 4 ) t , 0] + ... 

where the c's are determined from the d's by a rotational change 
of basis. Since homogeneous and isotropic interfaces do not 
display circular retardance or circular polarization, ct ( 3 ) is not 
included to simplify the mathematics. 

The C Vector for a Paraxial Ray 

The SAMEX instrumental polarization will be analyzed by a 
paraxial optics development. Consider a paraxial ray path 
through an optical system from surfaces q=l to Q with angles of 
incidence, i(q), and orientations of the plane of incidence, 

9(q). The Jones vector associated with the axial ray (down the 
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optical axis,) i(q)=0 for all q, is 

JJ = t (q) a(0) t(q-l) a( 0 ) . . . t ( 2 ) a(0) t(l) a(0) = T <r ( 0 ) 
where 

Q 

T - q I i t(q)- 

t 

The complex amplitude transmittance down the axis, T, is the 
product of the normal incidence complex amplitude transmittances 
at each surface. 

The Jones matrix associated with a ray at interface q can be 
expressed in terms of the expansion of the interface Jones matrix 
as 

JJ(i,v,q) = t(q) [c(0)+ i 2 (q ) (c( 0 , 2 ,q ) a(0 )+c( 1 , 2,q ) a( 1 )+c ( 2 , 2 ,q) a( 2) ) 

+ i 4 (q) (c( 0,4, q)a(0)+c( 1,4, q) a(l)+c( 2,4, q) a( 2) )+...) 

The Jones matrix associated with the entire paraxial ray 
path resulting from keeping terms to second order at each 
interface is (where X represents multiplication carried onto the 
next line) 

JJ = t(q)(o(0)+ i 2 (q)c(0,2,q)a(0)+c(l,2,q)a(l)+c(2,2,q)or(2)) X 

t ( q-1 )( a( 0 ) +i 2 ( q-1 ) (c(0, 2 ,q-l ) a( 0 )+c( 1, 2,q-l) o( l)+c( 2, 2,q-l) a( 2) )X. . 
t(2)( a(0)+i 2 ( 2) (c ( 0,2,2) a(0)+c( 1,2,2) o(l)+c( 2,2,2) a( 2) )X 
t(l)(a(0)+i 2 (l) (c(0, 2, 1) a(0)+c( 1,2,1) a( 1 )+c ( 2 , 2, 1 ) a(2) ) . 


Associated with each interface are four terms. Carrying out all 
the multiplications leads to 4^ terms, all in even powers in i. 
Collecting terms of equal power in i, there is one term at 
zero' th order and 3Q terms at second order. If i is assumed 
small, the large number of higher order terms are of diminishing 
importance. Collecting zero and second order terms in JJ Q and 
JJ 2 the expression for JJ is 

JJq + JJ 2 = T Cf(0) + T ff(0) I i 2 (q)c(0,2,q) 

q=i 

Q 2 

+ T a ( 1 ) Z i“(q)c(l,2,q) 

q=i 
Q 2 

+ T a(2) Z i 2 (q)c(2,2,q) . 

q=i 


Since no polarization or retardance was assumed on axis, the 
contributions to the second order polarization for this ray are 
j'ust sums of contributions from each surface. The multiplication 
taking place at second order for the a(l) term is of the form 


i 2 (q)<j(l) c(0) a( 0 ) a ( 0 ) . . . a( 0 ) + a(0) i 2 (q-l) <j( 1) a{0) a(0) . . . cr( 0 ) 

... + a( 0 ) u(0 ) <j( 0 ) a( 0 ) ...i 2 (q) 0 ( 1 ) 


where the c-dependence is not shown explicity. 

This equation contains the useful result that, when no 
elements display polarization or retardance at normal incidence, 
as in the SAMEX foreoptics, there is no order dependence in the 


second order terras. Only one non-identity matrix term occurs in 
each second order matrix product. The second order polarization 
associated with the paraxial ray path is obtained by a simple 
summation of second order polarization contributions at each 
intercept. Chipman (1987) gives a complete account of this 
derivation . 


Paraxial Optics Geometry 

The polarization aberrations for 5 AM EX are a description of 
the polarization behavior of an optical system expressed as an 
expansion about the center of the object and the center of the 
pupil. Thus it is appropriate and convenient to obtain the 
derivations from a paraxial ray trace; appropriate, because 
understanding the instrumental polarization near the center of 
the pupil and image is key to understanding instrumental 
polarization in general; convenient because the paraxial ray 
trace is linear, and thus easy to manipulate. 

The paraxial coordinate system used is a normalized right 
handed coordinate system. The z axis is the optical axis of a 
rotationally symmetric optical system. Light initially travels 
in the direction of increasing z. Figure C-l shows the notation. 

For a rotationally symmetric system, the object can be 
located on the y axis without loss of generality. The object 
coordinate H is normalized such that H = 0 in the center of the 
field (on the optical axis) and H = 1 at the nominal edge of the 


field of view. 

The location where a ray strikes the entrance pupil is 





PUPIL 



Figure C.l. Paraxial coordinate system. The paraxial system is 
a normalized right-handed coordinate system. The z axis is the 
optical axis of a rotationally symmetric optical system; -light 
initially travels in the direction of increasing z. Rays through 
an optical system are characterized by ray coordinates at the 
object and entrance pupil. H is the normalized object coordinate, 
p is the normalized pupil radius, and <fr is the polar angle in the 
pupil meas’Qlred counterclockwise from the y axis. The normalized 
Cartesian coordinates in the pupil are x and y. The chief and 
marginal rays are also shown. 
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specified by the polar pupil coordinates p and p is 
normalized such that at the edge of a circular pupil p 
The angle <p is defined here as it is in much of geometric optics, 
and in defiance to most analytical geometry, as being zero on the 
'y axis' and increasing counterclockwise. Normalized Cartesian 
pupil coordinates x and y can be used. They are defined as: 
x * p sin($) and y = p cos($). 

Expressions for the angle of incidence i and the orientation 

of the plane of incidence ^ of a ray at a given surface q will be 

expressed in terms of the marginal i(m,q) and _chief ray i(c,q) 

angles of incidence at that surface. Details of the derivation 

may be found in Chipman ( 1987, section D) . Note, however, for a 

radially symmetric system the angle of incidence should be a 
9 2 

function of H% p , and Hpcos<)» since the function should be . 
invariant to rotation of the system about the optical axis and 
must give the same result when -x is substituted for +x where 
x = pcos<j>. 

Assume that a paraxial ray trace has been performed for a 
specific system and that i(m,q) and i(c,q) have been 
calculated. A ray from normalized object coordinate H which 
passes through pupil coordinates p and <fi will have an angle of 
incidence i(q) and orientation of the plane of incidence 8(q) at 
surface q equal to: 

i(q) = [H 2 i 2 (c,q) + 2Hp cos((j>) i(c,q) i(m,q) + p 2 i 2 (m,q)] lj/ 2 
9 = sin -1 [p sin( $) i(m)/|i|] 


Figure C-2 shows the paraxial angle and plane of incidence 
for three field angles. The magnitude of the angle of incidence 
is represented by the length of the lines. The orientation of 
the lines corresponds to the orientation of the plane of 
incidence. Note that off axis, the pattern is a shifted version 
of the on axis pattern. The incidence angle is given by the 
pupil coordinate (p,$) and the image coordinate (H) since these 
coodinate define the optical path of a single ray. 


POLARIZATION ABERRATIONS 
Introduction 

This section derives the polarization aberrations for SAMEX 
as a Taylor series description of the instrumental polarization 
associated with paraxial rays through the foreoptics. 

Polarization aberration is a description of the polarization 
behavior of an optical system expressed by a expansion about the 
center of the object and center of the pupil. Table C-4 gives an 
overview of polarization aberration theory. The results are 
obtained in a form very similar to the wavefront aberrations. In 
particular, terms closely related to defocus, tilt, piston error 
as well as the Seidel and higher order aberrations can be 
associated with all eight of the basis Jones matrices. Since 
polarization effects are typically orders of magnitude smaller 
than wavefront effects, fewer terms are needed for a sufficient 
description. A method of calculating aberration coefficients for 
specific systems is developed in the next section. 
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Figure C.2. Paraxial angles and plane of incidence for three 
field angles. The magnitude of the angle of incidence is 
represented by the length of the lines. The orientation of the 
lines corresponds to the orientation of the plane- of incidence. 
Note that in the off-axis cases, the pattern is a shifted version 
of the on-axis pattern. > 




Table C-4. Polarization Aberration Theory Outline 


(1) Jones Vector 


The Jones vector is a complex Ex 

2-component vector describing J = ( . ) 

the electromagnetic field. Ey 

(2) Paraxial Ray Trace 

The angle of incidence, i, of any i(c), i(m) 

ray can be written in terms of the 
angle of incidence of the chief 
and marginal rays, i(c) and i(m). 

(3) Coating Calculation 

Thin film design programs t(s,i), t(p,i) 

calculate the amplitude 
transmission coefficients for s 
and p light, t(s,i) and t(p,i). 

(4) Taylor Series Representation 

The amplitude transmission 
coefficients, t(s,i) and t(p,i), 
are transformed to with a Taylor 
series expansion about normal 
incidence . 


t ( s , i ) a t Q (l + tfs,2)i 2 + ...) 
t (p, i ) = t Q (l + t(p, 2 ) i 2 + ...) 


(5) Jones Matrix for an Interface 


The second order Jones matrix for 
an interface can be written in terms 
of three Polarization Aberration 
terms represented by coefficients. 

(6) Aberrations Sum for System 

The second order polarization 
aberrations for the system is 
the sum of the aberration 
contributions of each interface. 


(7) Polarization Accuracy 
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P(l,2,0,0,q) 


P ( 1 , 0 , 2 , 2 ) = * 
P( 1 , 1,1,1) = g 
P( 1,2, 0,0) = g 


P( 1,0,2, 2,q) 
P(l,l,l,l,q) 
P( 1 , 2 , 0 , 0 ,q ) 


Summarize the performance of the system 
with a single number, the polarization 
accuracy (Ap). This is the maximum 
fraction of light coupled into the 
orthogonal polarization state. This 
occurs at the edge of the field of view. 
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The Polarization Aberration Expansion 

The wavefront polynomial expansion describes the variation of 
the optical path difference through an optical system as a function 
of ray coordinates. A closely related expansion will be presented 
for all four basis polarization matrices a(0), <j( 1), cr(2), a(3). 

The polarization aberration expansion for radially symmetric 
systems uses a very similar polynomial expansion to describe all 
eight basis polarization vectors. The principal difference is a 
modified form for the linear polarization and linear retardance 
terms since these involve both a magnitude and an orientation. 

The eight forms of polarization behavior can be 
characterized by four complex numbers, for example, the four 
elements of either the Jones matrix or the C vector. We 
introduce a new set of complex parameters, the complex 
polarization aberration coefficients which gives a description of 
the polarization behavior of an optical interface. It should be 
emphasized that the amplitude and phase of the coefficients are 
generally unrelated. They refer to different aspects of the 
instrumental polarization. The amplitude part of the coefficient 
describes amplitude and polarization effects while the phase part 
describes phase and retardance. When necessary the amplitude and 
phase of the polarization aberration coefficient, P, will be 
denoted by A and # where, 


P = A exp {7 $) . 


The complex aberration coefficient is written to contain 
polarization, amplitude, effects, retardance, and phase effects 
in a single term. The following polarization aberration 
expansion for zeroth and second order is used. The subscripts 
are defined as follows: 

P(k,u,v,w) = A(k ,u , v ,w) exp(T*(k,u,v,w) ) , 

where: k is the type of polarization behavior, 

u is the order of the H dependence, H u , 

v is the order of the p dependence, p v , and, 

w is the order of the $ dependence, (cos($)) w . 

The indices u, v, and w are used exactly as they are for the 

wavefront aberrations, as shown in the next section. 

The polarization aberration expansion of the Jones matrix 
for the SAMEX foreoptics is 

3 

JJ(H,p,<j>) = e c(k,H,,r,p) <j( k) 
k=0 

3 

- Z Z Z Z P(k,U,v,w) H u p v cos w (<J>) cr ( k ) . 
k=0 u v w 

Here we define the aberration coefficients as an expansion of the 
Jones matrix in terms of the ray coordinates p, $, and H. 

The C vector coefficients in this polarization aberration 
expansion are: 
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Amplitude and Phase term: 

c(0,H,p,*) = A(0, 0,0,0) + A(0,2,0,0)H 2 + 

+ A( 0,1, 1,1) H p cos ( $ ) + A( 0,0, 2,0) 

+ T( *(0,0, 0,0) + *(0,2,0,0)H 2 
+ *(0,1,1, 1)H p cos( *) + *(0,0, 2,0) p 

i 

Linear Polarization and Retardance terms: 

C(1,H,P,*) = A(l, 0,0,0) + A(l, 2,0,0) H 2 + 

+ H p (A(l, 1,1,1) cost*) - A(2, 1,1,1) 
+ p 2 (A(l, 0,2,0) cos ( 2 *) - A(2, 0,2,0) 

+ T (*(1,0, 0,0) + *(1,2, 0,0) H 2 + 

H p (*(1,1, 1,1) cos ( $ ) - *(2,1, 1,1) 
p 2 ( *(1,0, 2,0) cos(2*) - *(2,0, 2,0) 

Diagonal Polarization and Retardance terms: 

c(2,H,p,*) = A(2, 0,0,0) + A( 2 , 2 , 0 , 0 ) H 2 + 

+ H p (A(2, 1,1,1) cos ( * ) + A(l, 1,1,1) 
+ p 2 (A(2, 0,2,0) cos ( 2 i*) + A(l, 0,2,0) 
+ T (*(2,0, 0,0) + *(2,2, 0,0) H 2 + 

H p (*(2,1, 1,1) cos ( <j> ) + *(1,1, 1,1) 
p 2 (*(2,0, 2,0) cos ( 2 <*) - *(1,0, 2,0) 
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P 


sin ( *) ) 
sin (2*)) 

sin( *) )+ 
sin( 2 *) ) ) . 


sin( *) ) 
sin (2*)) 

sin( *) )+ 
sin( 2 *») )) 


Circular Polarization and Retardance terms: 





c(3,H,p,*) = A(3, 0,0,0) + A(3,2,0,0)H 2 + 

+ A(3, 1,1,1) H p cos(<j>) + A(3, 0,2,0) p 2 
+ T ( *(3,0, 0,0) + *( 3 , 2 , 0 , 0 ) H 2 
+ * (3,1,1,1)H p cos ( $ ) + 4( 3 , 0 , 2 , 0 ) p 2 ) , 

There are thirty two terms in this polarization aberration 
expansion to second order arising from four terms in each of the 
eight degrees of freedom of the Jones matrix. The terms may be 
grouped as follows: 


A( 0 ,u , v ,w ) 

Amplitude terms 

A( l,u,v,w) 

Linear polarization terms 

A( 2,u,v,w) 

Diagonal polarization terms 

A( 3 ,u , v , w) 

Circular polarization terms 

* ( 0 ,u , v ,w) 

Wavefront or phase terms 

4 ( l,u,v,w) 

Linear retardance terms 

* ( 2,u, v,w) 

Diagonal retardance terms 

* ( 3 ,u, v,w) 

Circular retardance terms 

P(k, 0,0,0) 

"Constant Piston" terms 

P(k, 2,0,0) 

"Quadratic Piston" terms 

P(k, 1,1,1) 

"Tilt" terms 

P(k, 0,2,0) 

"Scalar Defocus" terms 

P(k,0,2,2) 

"Vector Defocus" terms 


The names of the wavefront aberrations: piston, quadratic 


piston, defocus and tilt, are used here in an extended sense, to 


describe variations of components of the Jones vector which share 
the same functional dependences as the wavefront aberrations. 
Defocus is a p 2 variation of a parameter. Tilt is H p cos(*) 
variation. Quadratic piston is H variation. Thus, "amplitude 
defocus" means a p 2 amplitude variation. Likewise the "circular 
retardance tilt" is a the H p cos(*) circular retardance 
variation, and so on. 

This polarization aberration expansion is an equation which 
describes all possible second order variations of the Jones 
matrix, just as the second order wavefront aberration expansion 
spans the set of all second order wavefront variations. Thus the 
polarization aberration expansion characterizes quadratic 
variations of all forms of wavefront, amplitude, polarization and 
retardance . 

This polarization aberration expansion is a summation of 
terms in the different Pauli spin matrix components, not a 
product. Thus the four C vector elements can be pictured as 
acting in parallel, almost side by side in the aperture, but not 
in series. Each term describes an amount of a particular form of 
polarization, independent of the other contributions. 

An "aberration term" is to be considered as containing all 
the algebraic terms in the expansion with the same coefficient. 
Most of the coefficients occur only once and the aberration term 
contains only one algebraic term. The exceptions are the terms, 
A(l, 1,1,1), *(1,1, 1,1), A( 1 , 0 , 2 , 2 ) , *(1,0, 2, 2), A(2, 1,1,1), 

*(2,1, 1,1), A( 2 , 0 , 2 , 2 ) , and *(2,0, 2, 2). These aberration terms 
have components both along the axes and at 45 degrees . 


With SAMEX, the principal concerns are with the linear 
piston, linear tilt and linear defocus terms, both in 
polarization and retardance. These are going to be the largest 
terms present which corrupt the incident polarization state. 

These values are given in the Table 19 for SAMEX. 

For a detail discussion of the physical meaning of the 
polarization aberration coefficients see Chipman (1987), however 
a discussion of the orgin of tilt and piston, P(0,2,0), P(l,l,l), 
and P(2,0,Q) terms is included below. 

A distinction is made between scalar and vector 
aberrations. The wavefront aberrations are scalar aberrations, 
single valued functions of object and pupil coordinates. The 
linear polarization and linear retardance aberrations are vector 
aberrations since a magnitude and orientation is associated with 
these at each point. Amplitude, circular polarization and 
circular retardance aberrations are scalar since they are single 
valued and range positive and negative. 

Figure C-3 (top) shows the chief and limiting rays at an 
interface for objects on axis and at the edge of the field of 
view. Figure C-3 (bottom) is a plot of the value of the angle of 
incidence along the y axis as a function of p. Tilt terms 
naturally occur because as the object point moves off axis, the 
angle of incidence increases at one edge of the beam and 
decreases at the other edge. Tilt contains the first order 
portion of . this correction. 

Figure C-4 shows the off-axis angle of incidence squared and 
the decomposition of this into defocus, tilt and piston terms. 


ON-AXIS 



ANGLE OF li 
INCIDENCE 



B. 


Figure £3 Angles of incidence for objects on and off axis. I 
the top frame, the chief and limiting rays at an interface are 
shown for objects on the optical axis and at the edge of the 
field of view. In the lower part of the figure, the angle of 
incidence for rays incident along the y axis (in the paraxial 
system) is plotted as a function of the normalized radius p. 
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Figure C.Jp Quadratic effects of off-axis angles of incidence. 

The squafe of the off-axis angle of incidence is shown along with 
its decomposition into defocus, tilt and piston terms. 



These terms are required to describe a quadratic variation whose 
vertex is located at an arbitrary position on the y axis because 

x 2 + (y-a) 2 = x 2 + y 2 - 2y a + a 2 . 

In this case since a is a linear function of H, 

a * k H , 

then the quadratic polarization variation becomes 

c(x 2 + (y-kH) 2 ) = c(x 2 +y 2 ) - 2c y k H + c k 2 H 2 

■ P(0,2,0)p 2 +P(1,1,1)2H p cosU) + P(2,0,0)H 

where P(0,2,0), P(l,l,l), and P{2,0,0) are the defocus, tilt and 

quadratic piston aberration coefficients and the polarization 

index is not shown. Tilt and piston terms arise naturally from 

decentered defocus. Similarly, the fourth order wavefront 

aberrations coma, astigmatism, field curvature and distortion 
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arise naturally from decentered spherical aberration, p . 


CALCULATION OF ABERRATION COEFFICIENTS 

The method used to calculate the second order polarization 
aberration coefficients for the SAMEX foreoptics given the C 
vector power series for each interface is detailed in this 
section . 

Single Surface Aberrations for Amplitude Transmittance Relations 
For homogeneous and isotropic interfaces characterized by 
amplitude transmittance relations , such as lenses, mirrors and 
thin film coatings, the polarization aberrations at a interface 
simplify considerably. At these interfaces the Fresnel equations 
and related thin film equations are separable into s and p 
components, so the Jones matrices representing the interface in 
s-p coordinates are diagonal. The off-diagonal terms, diagonal 
polarization cr(2) and circular polarization a(3), are not 
present. Further, with isotropic media, the s and p amplitude 
transmission coefficients at normal incidence must be equal. 

Thus the amplitude transmission functions can be expanded as: 

a(s,i) = a(0)(l + a(s,2)i 2 +...) exp(T(d(0) + d(s,2)i 2 +...) 

* t ( l+(a(s, ?) + T d(S,2))i 2 +...) , 

a(p,i) = a ( 0 ) ( 1 + a(p,2)i 2 +...) exp(T(d(0) + d(2,p)i 2 +...) 

= t (1 + a(p,2) + T d(p,2)i 2 +...) 


where : 


t 


a(0) exp( 7 d(0) ) 


The s-p coordinate Jones matrix expansion to second order is 


JJ(i) 


a ( s , i ) 



1^0 a ( p, i ) 


= t ( cr( 0 ) ( 1 + A( 2 ) + Td( 2 ) ) i 2 ) + o(l)(a(2) + Td ( 2 ) ) i 2 > . 


where 


A( 2) 

= (a(s,2)+a(p,2) )/2 

a( 2) 

= (a(s,2)-a(p,2) )/2 

D(2) 

= (d(s,2)+d(p,2) )/2 

d ( 2 ) 

= (d(s,2)-d(p,2) )/2 


The s-p coordinate C vector Taylor series expansion to second 
order is 


C ( i ) = t (1,0, 0,0) + i 2 t (A( 2)+TD( 2) , a(2)+Td(2), 0, 0). 


The x-y coordinate C vector Taylor series for orientation of the 
plane of incidence 0 is 


h; 


C(i,8) = t(l, 0,0,0) + i 2 t (A( 2 ) +~D( 2 ) , 

( a ( 2 ) +7d ( 2 ) ) cos (20), (a ( 2 )+id ( 2 ) ) sin ( 29 ) , 0). 

The normal-incidence polarization aberration terms (the constant 
piston terms) are zero: 

P(l, 0,0,0) = P(2, 0,0,0) = P(3, 0,0,0) = 0 . 

There is no polarization or retardance on axis, only the 
amplitude and phase transmission factor t. 

All terms for off-axis diagonal and circular polarization 
are zero: 

P(2,u,v,w) = P(3,u,v,w) = 0 

Thus, the single surface C vector in paraxial coordinates is 
obtained by substituting the paraxial representation of i(H,p,$) 
and v(H,p,<j>) into C(i,v) yielding 

c(0,H,p,+) = t + t d(0,2) [ H 2 i 2 (c) 

+ 2H p cos($) i(c)i(m) + p 2 i 2 (m)] 

c(l,H,p,$) = ' t d( 1,2) [ H 2 i 2 ( c ) 

+ 2H p cos(<^) i(c)i(m) + p 2 cos(2^) i 2 (m)] 

c(2,H,p,4) = t d(l,2) [ 2H p sin(^) i(c)i(m) 

+ p 2 sin(2^) i 2 (m)] 


c(3,H,p,<fr) = 0. 


Since there is no diagonal polarization, the only contributions 
to c(2) arises from the rotation of linear polarization from the 
s-p coordinates into the x-y coordinates. 

» 

Polarization Aberration Coefficients for Systems 

Since the polarization aberrations are only being evaluated 
to second order in the angle of incidence, the difference between 
spheres, parabolas, conics or other radially symmetric aspherics 
does not occur at this order. The relevant shape parameter here 
is only the vertex radius of curvature. The angle and plane of 
incidence differences between these types of interfaces are the 
same at second order but will differ at fourth order and higher. 

For surfaces q=l to 0, each surface is characterized by 
three complex parameters from the normalized C vector expansion: 


. d(0,0) = t(q), d ( 0 , 2 ,q) , and d(l,2,q) . 


The single surface polarization aberration coefficients are: 


P(0,0,0,0,q) 
P(0,2,0,0,q) 
P(0,l,l,l,q) 
P(0,0,2,0,q) 
P(l,2,0,0,q) 
P( 1, 1,1,1, q) 
P(l,0,2,2,q) 


t(q) 

t ( q ) d( 0 , 2 ,q ) i 2 (c) 

2 1 ( q ) d( 0 , 2 ,q) i(c)i(m) 
t(q) d(0,2,q) i 2 (m) 
t (q ) d( 1 , 2 ,q) <i 2 (c) 
2t(q) d(l,2,q) i(c)i(ra) 
t (q) d( 1 , 2 ,q) i 2 (m). 


The polarization aberration coefficients for the system are 
calculated by chain multiplying the single surface polarization 
aberration expressions and keeping terms to second order in H and 
p. The zero and second order Jones matrices for the q'th 
interface are: 


= P(0f0f0,0,q) a(0) = d(0,0,q) cr(Q) = t(q) a ( 0 ) , 

JJ 2 (q,H,p,<») = a(0) [H 2 P(0,2,0,0,q) + 2H p cos ( *) P ( 0 , 1 , 1 , 1 ,q ) 

+ p 2 P(0,0,2,0,q) ] 

+ a ( 1 ) [ H 2 P(l,2,0,0,q) + 2H p cos ( * )P ( 1 , 1 , 1 , 1 ,q ) 

+ p cos ( 2 ) P { 1 , u , 2 , 2 , q ) ] 

+ a ( 2 ) [ 2H p sin(*)P(l,l,l,i,q) + p 2 sin( 2*)P( 1, 0 , 2 , 2,q) 1 
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= a ( 0 ) t(q)d(0,2) [ H 2 i 2 (c) + 2H p cos ( 4 ) i ( c ) i (m ) 

+ p 2 i (m) 2 .] 

+ a(l) t(q)d(l,2) [H 2 i(c) 2 + 2H p cos( <p) i (c) i (m) + p 2 cos(2$) 

t 

+ a(2) t(q)d(l,2)[ 2H p sin( <p) i (c ) i (m) + p 2 sin(2$) i 2 (m) ] 

Multiplication of the single surface Jones matrices yields 
JJ(H,r,p) = 11 JJ(q,H f p,<p) 


q-Q,-l 

q=Q,-l. C JJ 0 (< 3 } + ] * 


Since JJ Q (q) is a constant function, independent of H, p, and <p , 
the { H , p , $ ) dependence can be dropped. This expression contains 
20 terms. The order of a term is the sum of the powers of H 
and p , H u p v cos(<p) w , i.e. order = u + v 


There is one first order term and Q second order terms. 


i 2 ( 


The zero order Jones matrix is 



H 


JJ 0 = q=Q,-l ^iq) - T / 

the system amplitude transmittance. The second order Jones 
matrix is greatly simplified since, for isotropic surfaces, all 
zeroth order Jones matrices are a constant times the identity 
matrix <j( 0). The second order only contains products which 
contain a single second order term.’ The second order Jones 
matrix is 

JJ2(H,p,$) = JJ2 ( q, H , p , $ ) 

At second order the weakly polarizing isotropic interfaces do not 
display order dependence. The product of any two second order 
terms is fourth order. The order dependence enters at fourth and 
higher order. Second order is a simple sum of polarization 
contributions. Collecting the piston, tilt and defocus terms 
from the second order Jones matrix provides the coefficients for 
the system polarization aberration expansion to second order: 


P(0, 0,0,0) 

S 

T 

t 



P(0,2,0,0) 

= 

T 

z 

q=l,Q 

d{0,2,q) 

i 2 (c) , 

P(0, 1,1,1) 


2T 

Z 

q=l,Q 

d(0,2,q) 

i(c) i(m) 

P(0, 0,2,0) 

= 

T 

Z 

q=l,Q 

d(0,2,q) 

i 2 (m) , 

P(l,2,0,0) 

= 

T 

1 u t-i 
O 

d( 1 , 2 ,q) 

i 2 (c) , 

P(l, 1,1,1) 


2T 

L 

q=l,Q 

d{ 1 , 2 ,q) 

i (c) i (m) 

P( 1 , 0 , 2 , 2 ) 

=: 

T 

Z 

q=l,Q 

d(l,2,q) 

i 2 (m) 



The other three zero order coefficients and the other six second 
order coefficients are all zero: 


P(l, 0,0,0) 
P( 2 , 2 , 0 , 0 ) 
P(3,2,0,0) 


P(2, 0,0,0) 
P(2, 1,1,1) 
P(3, 1,1,1) 


P(3, 0,0,0) = 0 , 

P(2, 0,2,0) = 0 , 

P(3, 0,2,0) = 0 . 


The amplitude and polarization coefficients are the real 
parts of the P coefficients 


A(k,u,v,w) * Re(P(k,u,v,w) ) 


The retardation coefficients are the imaginary parts 


*(k,u,v,w) = Im(P(k,u,v,w) ) 


The polarization aberration coefficients are calculated for the 
foreoptics from the paraxial geometry and the normalized C- 
vectors for the coatings and interfaces. 

Polarization Accuracy 

The Jones vector gives the amplitude of the electric field 
and the square of the amplitude gives the intensity of the 
components. With the Jones matrix one is able to calculate the 
polarization effects for the optical system. From above the 
second order the Jones matrix which gives the linear polarization 
and linear retardence is 


JJ 9 (H,p,*) = <x(l)[H* P(l,2,0,0) + 2 H p cos(*) P(l, 1,1,1) 


+ p* cos ( 2 $) P( 1 , 0 , 2 , 2 ) ] 


+ ff(2)[(2Hp sin(*) P(l, 1,1,1) + p 2 sin(2*)P(l,0,2,2)] 


The on axis linear polarization and linear retardance of the 
SAMEX foroptics are characterized by the term linear defocus. 

The instrumental polarization function JJ(H,p,<») for linear 
de focus is 


JJ(H,p,$) = t[<Jq + P(l,Q,2,2)p (<jjCOs2ifr - C2sin2<fr] 


= t [ a n + (a 


1022 


+ j 6 ^q 22 ) ( <I ^cos2^ - a 2 s *- n 2$)l 


Here, r is the amplitude transmittance of the system down the 
optical axis. It (t) describes the polarization independent 
reflection and absorption losses associated with the ray down the 
optical axis at normal incidence at all interfaces. P(l,0,2,2) 
describes the linear polarization (^1022^ and linear retardance 
(<Siq 2 2 ) associated with the marginal ray. 

Maximum coupling occurs when the incident light is 
circularly polarized, since circularly polarized light can always 

A A 

be decomposed into equal components of J r and J fc everywhere in 
the pupil. The coupling is zero in the center of the pupil 
(where the polarization and retardance vanishes) and increases to 
a maximum coupling fraction of 


I (H,l,$) = 

c f max 


P(l,0,2,2) | 2 = a * , + 5 


1022 


1022 


at the edge of the pupil. The net fraction of incident circular 
polarized light coupled into the orthogonal circular polarized 
state is given by the integral over the pupil of 


2 2n i 

x c = -L—L- / d* / pdplp tP(1,0,2,2) 1 2 

* 0 0 ' 

_ | tP( 1,0,2,2) l 2 
3 


For incident linear or illiptical polarized light, the fraction 
of coupled intensity is less because the light is not composed of 
equal fractions of eigenstates. 

The coupling is minimum for incident linear polarized light, 
which will be in one of the eigenpolarizations along one axis in 
the pupil and in the orthogonal eigenpolarization along the 
orthogonal axis. The fraction of coupled energy will be 
calculated assuming an incdient polarization state of horizontal 

A 

linear polarized light H for calculational simplicity, the same 
fraction is coupled for any linear polarized incident state. The 
polarization state transmitted by an optical system described by 

A 

linear polarization defocus for H in is 



The fraction of incident H light couplied into V light is equal 
to 


I 


c 



2w 

/ d<j» 
0 



2 TT 

/ d* 
0 


= | tP(1,0,2,2) | 2 

TT 


. | tP( 1 # 0 y 2 # 2) | 2 
6 


/ pdp|3<H,p,$) . V 


/ pdp|P(l,0,2,2)p I 2 sin2$| 2 


2 IT • 2 , 1 e 

/ sin 2$d$' / p ap 

0 0 


This is the mininmum fraction of energy coupled by linear defocus 
aberrations. Since any elliptically polarized incident beam can 
be written as a sum of linear and circularly polarized light, the 
coupling fraction for arbitrarily polarized light lies in the 
range 


I T P( 1 , 0 , 2 , 2 ) I 2 

6 


xP ( 1 , 0 , 2 , 2 ) 
3 


When unpolarized of circularly polarized light is incident 
on the optical system, maximum polarization coupling occurs. 
This maximum is the polarization accurracy of a system and is 
calculated by the coupling integral^ ^cmax' 


I 

c ,max 



pupil 


c 2 ( 1 ) + c 2 (2) + c 2 ( 3 ) | pdpd*. 



This integral must be evaluated numerically except for some 
special cases (see section 3.5). However, an analytic upperbound 
on the polarization accurracy can be easily established by using 
the triangle inequality 


c,raax ‘ 1 I JpupU I((Re[c(1,1)2 + <Re[c(2)] 2 + (RelcU]) 2 ) 2 
+ < (Im[c(l)] ) 2 + ( Im [c( 2 ) ] ) 2 + ( Im [c ( 3 ) ] ) 2 ) 2 ] pdpd* . 


We can evaluate the upper bound on the polarization accurracy 
using the second order polarization aberration expansion 
coefficients. The following integration is for the real part of 
the aberration coefficients. a luvw 


I T | 2 2ir 1 2 22 

I 58 / d<fr / pdp[(a 120Q H + a im H P cos, fr + a i022 p cos2 *> 


2 2 

+ (a Hll Hpsin, < ) + a 1022 p sin2l * ,) 1 


2 2 2 

I 1 2 r _ 2 „4 , a llll H . a 1022 , 

“ I T I a i200 H + 2 + T~ ^ * 


The corresponding integral, imaginary ' * s obtained by 
substituting the imaginary part of the polarization aberration 
coefficients, 5 luvw for a luvw in the above expression. The 
upperbound on the polarzation coupling is 



I c / max < ^real + imaginary 


2 2 2 2 2 2 
2 r 2 „4 a llll H a 1022 P 2 „4 . 5 1111 H (S 1022 1 

T| l a i?nn H + 5 + — T — + ° 1200 H 3 + — T — J 


- | t| 2 [ |p< i, 2 ,o,o) | 2 h 4 . .|m,i,i.i)|V „ ]_ 


This is the second order upper bound to the polarization coupling 
for systems of weak polarizers. 

Hence the average effect over the image and pupil can be 
obtain by integrating over p, and $. The polarization 
accuracy, A^, defined as the maximum fraction of light 
(intensity) which can be coupled into orthogonal polarization 
state. The incident polarized state is given by the Jones 
vector, 3. This vector is effectively rotated by the optical 
system and the rotation is given by the Jones matrix, JJ. The 
amount of polarization along the orthogonal state of 


polarization, of the incident polarization state is g 

the projection of JJ(5) into J', 


- . .»a r> Ki r 

^ x v u ujr 


JJ( 3) -5' . 


This value is given by the square of the second order Jones 
matrix and is given in terms of the polarization aberration 
coefficients. (For the incident light in the polarization 

A 

state I Q p, then the result of instrumental polarization is to 
couple the polarization into the orthogonal state q. The amount 


of coupling is given by 


Ip I_(l-A )p + Alq 
o o p p o 

where A p = polarization accuracy.) From the above results, the 
polarization accuracy is given in terms of the second order 
polarization abberation coefficients P(l,0,2,2), P(, 1,1, 1,1), and 
P( 1,2, 0,0) : 


A p = j / |jJ 2 (H = l,p,«M| 2 pdpd* 
or 

A p = 2P 2 (1,2,0,0) + -|p 2 (1, 1,1,1) + -jP 2 ( 1 , 0 , 2 , 2 ) 

where the integrations have been carried out for the squared 
terms and estimated for the crossed terms. For the SAMEX 
magnetograph design given herein the polarization accuracy value 
of 


A < 1.4 x 10 ' 

P 

is obtained for specially selected optical coatings was 

obtained. These second order aberration coefficients are given 

in Section II-4. The second order coefficients are sufficients 

since the next order that contributes is the forth order. The 

2 

polarization effects would be on the order of ( A ) . 

P 
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